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ABSTRACT 


J 

An  analytical  and  experimental  investigation  was  made  of  the  aeroelastic 
flutter  and  divergence  behavior  of  a  forward  swept  graphite/epoxy  wing 
aircraft  with  rigid  body  freedoms  in  free  flight. 

Analytically,  a  transient  aeroelastic  analysis,  using  a  Rayleigh-Ri tz 
formulation,  modified  strip  theory  and  measured  aircraft  aerodynamic 
derivatives,  was  developed  to  predict  and  interpret  the  aeroelastic  behavior 
and  instability  mechanisms.  The  transient  aerodynamics  approximated  the 
Theodoresen  function  by  Pade  approximants  and  used  augmented  state  variables 
to  investigate  stability  by  standard  linear  eigenvalue  analysis.  Some  effects 
of  aerelastic  tailoring,  wing  frequencies,  rigid  body  modes,  static  stability 
margins,  and  inertia  coupling  effects  were  explored. 

Experimentally,  a  generic,  full-span,  30f  forward  swept  wing  aircraft 
model  allowing  both  rigid  body  pitch  and  plunge  was  constucted  and  tested  in 
MIT's  low  speed  wind  tunnel.  The  effects  of  aeroelastic  tailoring  were 
demonstrated  by  using  wings  with  different  ply  layups-^  namely  [0-/90]  , 
[15-/0]  ,  [30-/0]  ,  and  [-15-/0]  y?which  covered  a  range'of  bending-twisting 
couplinf.  The  win§  tunnel  tests  ofi  these  softly  restrained  .^free  flying"' 
models  Revealed  body  freedom  flutter,  bending-torsion  wing  flutter,  and  a 
tunnel  support  related  dynamic  instabilty  which  could  be  eliminated  by  proper 
adjustment  of  the  support  stiffness.  Good  agreement  with  linear  theory  was 
found  for  all  observed  instabilities  and  transient  responses.  The  body 
freedom  flutter  boundary  followed  the  same  aeroelastic  tailoring  trends  as  the 
cantilever  divergence  boundary  but  at  a  lower  speed.  Additional  tests  with 
only  rigid  body  pitch  present  showed  still  lower  flutter  speeds.  The  model 
with  [15-/0]  wing  showed  the  best  tailored  aeroelastic  behavior  of  all  wings 
tested.  S 
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FOREVARD 

This  report  describes  work  done  at  the  Technology  Laboratory  for  Advanced 
Composites  (TELAC)  of  the  Massachusetts  Institute  of  Technology  for  the  Air 
Force  Office  of  Scientific  Research  under  Contract  No.  AFOSR 
F49620-85-C-0099.  Dr.  Anthony  K.  Amos  was  the  technical  monitor. 

The  work  reported  here  was  performed  during  the  period  1  July  1985 
through  30  June  1986,  and  represents  an  Sc.D.  thesis  by  Gun-Shing  Chen 
entitled,  "Aeroelastic  Behavior  of  Forward  Swept  Graphite/Epoxy  Uing  Aircraft 
with  Rigid  Body  Freedoms",  May  1986.  The  work  was  done  under  the  supervision 
of  John  Dugundji,  the  Principal  Investigator,  and  included  the  assistance  of 
an  undergraduate  student,  Karen  Needels,  and  the  supporting  laboratory  staff. 

A  brief  paper  based  on  some  of  this  work  entitled,  "Experimental 
Aeroelastic  Behavior  of  Forward  Swept  Graphite/Epoxy  Wings  with  Rigid  Body 
Freedoms",  was  presented  at  the  27th  AIAA/ASME/ASCE/AHS  Structures,  Structural 
Dynamics  and  Materials  Conference,  San  Antonio,  Texas,  May  19-21,  1986,  AIAA 
Paper  No.  86-0971,  and  will  be  published  shortly  as  an  article  in  the  Journal 


of  Aircraft. 
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CHAPTER  ONE 


INTRODUCTION 

In  the  early  history  of  forward  swept  wing  (FSW)  design, 
e.g.,  Ju-287  bomber,  HFB-320  corporate  jet  and  XFG-1  fuel 
transport  glider,  the  wing  was  swept  forward  moderately  for 
non-performance  related  purposes  and  no  serious  aeroelastic 
problem  arose  within  their  designated  flight  envelope.  For  the 
high  speed  performance  FSW  aircraft,  the  adverse  aeroelastic 
behavior,  presumed  aeroelastic  wing  divergence,  became  critical 
and  virtually  impossible  to  improve  without  incurring  a 
significant  structural  weight  penalty.  Since  then,  this 
seemingly  insurmountable  aeroelastic  problem  excluded  the  FSW 
concept  from  practical  aircraft  design.  It  was  not  until  the 
1970's  that  the  FSW  concept  was  revived  by  the  successful 
development  and  application  of  advanced  composite  materials  to 
conventional  aircraft.  Its  application  to  the  aeroelastic 
problem  of  FSW  design  was  first  studied  by  Krone  in  the  mid 
1970's  (1).  His  work  turned  the  key  to  unlock  the  FSW's 
inherent  aeroelastic  problem.  Since  then,  renewed  interests  in 
FSW  design  have  inspired  many  studies  and  tests  concerning  the 
aeroelastic  behavior,  aeroelastic  tailoring  [2-12],  and  active 
aeroelastic  instability  control  technology  [13,14]  for  such 
wing  configuration.  The  X-29  experimental  aircraft 
demonstrates  the  most  visible  revisit  of  the  FSW  design 
concept . 


Because  of  the  inherent  aeroelastic  destiffening 
characteristics  of  FSW  configuration,  it  is  natural  to  assume 
that  this  design  would  be  divergence  critical.  However,  when 
the  fuselage  is  allowed  to  participate  in  the  motion,  the 
situation  can  be  markedly  different.  An  aeroelastic/f light 
dynamic  interaction  becomes  strongly  possible  due  to  the 
attendant  aerodynamic  stiffening  phenomenon  of  the  aircraft 
short  period  mode.  Such  interaction  may  modify  the  presumed 
cantilever  wing  aeroelastic  behavior  or  even  cause  additional 
aeroelastic  problems.  Therefore,  the  effect  of  rigid  body 
freedoms  on  the  aeroelastic  behavior  of  a  FSW  aircraft  is  more 
critical  than  it  is  for  conventional  aircraft. 

One  of  the  earliest  reports  on  the  subject  of  rigid  body 
freedom  effects  was  due  to  Frazer  and  Duncan  in  1929  [15]. 
Their  study  showed  that  fuselage  mobility  had  a  negligible 
effect  on  the  symmetric  wing  flutter  speeds,  but  significantly 
increased  the  antisymmetric  wing  flutter  speed.  No  flight 
dynamics  involved  instability  was  mentioned  in  their  report. 

In  the  mid  1950s,  Gaukroger  [16,17]  conducted  a  series  of 
wind  tunnel  investigations  concerning  the  effect  of  wing  root 
boundary  conditions.  His  investigations  showed  that  the  rigid 
body  effect  exists  only  for  the  "flying  wing"  type  aircraft 
with  a  low  value  of  the  fuselage  pitching  moment  of  inertia, 
but  rarely  exists  for  a  pratical,  conventional  aircraft.  These 
investigations  of  aft  swept  and  delta  wing  configurations  were 
summarized  in  Ref.  18.  The  term  "body  freedom  flutter"  was 
then  commonly  used  to  characterize  the  instability  resulting 


from  wing  bending  coupling  with  rigid  body  pitching  and 
plunging  freedoms,  it  was  also  mentioned  that  symmetric  wing 
flutter  is  barely  influenced  by  the  presence  of  rigid  body 
freedoms . 

In  1954/  McLaughlin  performed  a  theoretical  investigation 
of  short  period  dynamic  instability  with  an  aft  swept  elastic 
wing  [19].  He  found  no  dynamic  instability  due  to  wing 
flexibility.  In  1955,  Cunningham  and  Lundstrom  also  described 
the  flutter  testing  of  a  rocket  mounted  straight  wing  [20]. 
Their  observations  revealed  a  low  frequency  type  instability, 
involving  rigid  body  motions  and  wing  deformation,  in  addition 
to  the  wing  flutter.  For  the  antisymmetric  oblique  wing 
aircraft  [21],  Jones  and  Nisbet  had  demonstrated  that  rigid 
body  freedoms,  including  roll  freedom,  are  essential  to 
aeroelastic  analysis.  With  the  proper  ratio  of  wing  roll 
inertia  to  fuselage  inertia,  a  low  frequency  body  freedom 
flutter  occurred  at  a  speed  greater  than  the  wing  divergence 
speed . 

Due  to  the  presumed  static  aeroelastic  wing  divergence, 
the  aircraft  rigid  body  freedoms  were  usually  ignored  in  many 
of  the  renewed  FSW  studies  prior  to  1980.  Recently,  however, 
several  authors  pointed  out  significant  effects  of  rigid  body 
freedoms  in  modifying  the  cantilever  aeroelastic  behavior  of 
FSW  design  .  While  performing  studies  for  a  Rockwell  design  of 
an  FSW  demonstrator,  Miller  and  Wykes  discovered  a  low 
frequency  type  instability  involving  a  coupling  between  the 
wing  divergence  mode  and  the  aircraft  short  period  mode  [22]. 
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Their  analysis  showed  that  the  rigid  body/wing  bending  flutter 
depends  on  the  distance  between  the  wing  center  of  mass  and  the 
aircraft  center  of  mass  for  a  typical  FSW  aircraft.  In  Ref.  23 
this  effect  was  explored  on  a  limited  basis  by  a  simple 
half-plane  model  without  rigid  body  plunging  motion  present. 
Further  work  on  a  large,  half-plane  aircraft  model  was 
performed  by  Chipman  et  al  [24,25]  at  NASA  Langley.  Generally, 
however,  the  aeroelastic  behavior  of  a  free  flying  FSW  aircraft 
is  not  yet  well  understood,  and  experimental  data  on  this 
effect  is  limited. 

The  mechanism  by  which  cantilever  wing  divergence  could  be 
eliminated  or  delayed  through  the  use  of  advanced  composite 
materials  is  based  on  the  inherent  property  of  advanced 
composite  materials,  namely  directional  stiffness.  The 
application  of  directional  stiffness  in  structural  design  for 
deformation  control  under  aerodynamic  loading,  or  aeroelastic 
tailoring,  is  not  totally  new  in  aeronautical  history.  This 
concept  was  applied  as  early  as  1949  by  Munk  in  a  wooden 
propeller  design  (26).  The  intended  fixed  pitch  in  Munk's 
design  was  accomplished  by  properly  orienting  the  grain  of  the 
wood.  The  other  early  example  was  a  novel  wing  design,  known 
as  the  aero-isoclinic  wing,  incorporated  in  the  design  of  Short 
Sherpa  prototype  [27].  Aeroelstic  tailoring  in  this  case  was 
achieved  in  part,  by  placing  the  torsion-box  well  back  in  the 


wing 


For  the  past  three  decades,  however,  the  concept  of 


aeroelastic  tailoring  was  not  fully  realizable  until  the 


development  of  advanced  composite  materials,  which  have  the 
advantage  of  easily  orienting  the  stiffness  properties,  as  well 
as  high  strength  and  stiffness  to  weight  ratios.  Following  the 
successful  application  of  advanced  composite  materials  to  the 
primary  structure  of  conventional  aircraft,  such  as  AV-8B  and 
F-18,  the  concept  of  aeroelastic  tailoring  using  advanced 
composite  materials  was  first  applied  to  a  remotely  piloted 
research  vehicle  [28],  HiMAT,  in  1978.  As  in  the  previous 
example,  the  purpose  of  the  aeroelastic  tailoring  in  HiMAT  is 
to  maintain  the  desired  aerodynamic  shape  and  aerodynamic 
performance  at  varied  flight  regimes.  It  is  only  in  the  case 
of  a  forward  swept  wing,  that  aeroelastic  tailoring  is  not  only 
designed  for  aerodynamic  shape  control  but  also  aiming  at  a 
control  of  aeroelastic  instability 

In  previous  investigations  at  MIT,  the  aeroelastic 
behavior  of  unswept  and  forward  swept  graphite/epoxy  cantilever 
wings  was  studied  experimentally  [9,10].  The  present 
investigation  extends  these  investigations  to  include  the 
effect  of  rigid  body  freedoms  on  the  aeroelastic  behavior  of 
FSW  aircraft.  The  principal  objectives  in  this  investigation 
are  to  study  the  aeroelastic  behavior  of  FSW  aircraft  and 
aeroelastic  tailoring  with  advanced  composite  materials.  The 
research  effort  consists  of  extensive  "free  flying"  wind  tunnel 
experiments  and  analytical  investigation.  A  linear  aeroelastic 
analysis  using  modified  strip  theory  and  measured  aircraft 
aerodynamic  derivatives  is  developed  to  predict  and  interpret 
the  expremental  aeroelastic  behavior.  The  aeroelastic  behavior 


of  the  FSW  aircraft  is  then  carefully  examined  through 
parametric  studies  of  the  verified  analytical  model.  The 
graphite/epoxy  wings,  with  various  degree  of  bending-twisting 
coupling,  are  formulated  by  a  Rayleigh-Ri tz  approximation.  The 
transient  aeroelastic  analysis  is  accomplished  by  a  rational 
polymoninal  fit  of  the  Theodorsen  function  in  the  Laplace 
domain.  Stability  of  the  wind  tunnel  support  system  is 
carefully  analyzed  for  the  "rigid"  model.  Descriptions  of  the 
af ormentioned  analytical  formulation  are  contained  in  Chapter 
2. 

In  Chapter  3,  details  are  given  of  the  experimental 
design,  setup  and  procedure  used  to  detect  the  aeroelastic 
behavior  of  the  FSW  aircraft  model.  The  results  of  the 
experiments  are  given  in  Chapter  4.  A  discussion  concerning 
the  correlation  of  observed  aeroelastic  behavior  and  predicted 
results  are  included  in  Chapter  5.  Also  in  Chapter  5,  the 
criticality  of  body  freedom  flutter  instability  is  compared 
with  the  cantilever  wing  divergence  instability.  The  effect  of 
aeroelastic  tailoring  on  body  freedom  flutter  is  examined  for 
all  wings  tested.  The  mutual  effect  of  support  instability  and 
body  freedom  flutter  is  examined  analytically  and 
experimentally.  Lastly,  parametric  trends  are  predicted  for 
aircraft  static  margin,  ratio  of  wing  pitch  inertia  to  fuselage 
inertia,  and  wing  frequency.  Chapter  6  summarizes  the 
conclusions  and  recommendations  of  this  investigation. 
Appendices  are  given  which  detail  the  analytical  formulation 
and  wind  tunnel  experimental  results. 
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CHAPTER  TWO 
THEORY 

2 . 1  Overview 

In  the  present  study,  the  aeroelastic  tailoring  concerns 
the  application  of  advanced  composite  materials  to  the  control 
of  aeroelastic  behavior.  From  the  underlying  classical 
laminated  plate  theory  (29],  the  out  of  plane  behavior  of  a 
symmetric  laminate  can  be  analytically  described  by  the 
laminate  constitutive  relation  as 

Du  Dn  Di6  ]  (  *x 

°12  °22  °26  Ky 

°16  d26  °66  J  1  Kxy 

For  the  moderate  to  high  aspect  ratio  wing,  the  chordwise 

rigidity  is  usually  assumed.  In  this  case,  the  most  important 
parameter  in  the  aeroelastic  tailoring  analysis  is  the  D^g  term 
which  dictates  the  bending- twi sting  coupling,  i.e.,  Mx  to  Kxy 
and  Mxy  to  «x  behavior,  of  a  laminated  wing  structure.  Usually 
this  D^g  term  is  normalized  by  the  bending  stiffness,  D^,  and 
torsional  stiffness,  Dgg,  as  Dig//^I7  Dgg  to  indicate  the 
effect  of  elastic  coupling  on  aeroelastic  tailoring  (8). 
Physically,  this  coupling  stiffness  ratio  is  bounded  in  between 
-1  and  1  which  can  be  easily  proved  from  the  positiveness 


condition  of  the  bending  stiffness  matrix  in  E q.  (2.1).  For 
wings  with  plate-like  deformation,  e.g.,  low  aspect  ratio  or 
highly  elastically  coupled  wing,  it  should  be  noted  that  the 
D2g  term  will  play  an  increasingly  important  role  in  the 
aeroelastic  tailoring  mechanism. 

A  simple  example  of  elastic  bending-twisting  coupling  is 
given  in  Figure  2.1.  The  degree  of  tailoring,  or 
bending-twisting  coupling,  is  considered  to  be  dominated  by  the 
designated  major  ply  orientation  as  shown  in  Figure  2.1.  In 
comparing  with  the  wing  structural  axis,  (laminate  axis). 
Figure  2.1. a  has  the  major  ply  orientation  coinciding  with  the 
laminate  axis,  and  leads  to  a  vanishing  or  negligible  D^g 
term.  Figure  2.1.b  and  2.1.c  have  the  major  ply  orientation 
rotate  toward  and  away  from  the  wing  leading  edge  respectively, 
and  will  consequently  have  favorable  and  adverse  D^g 
contributions.  The  deformation  coupling  will  be  illustrated  by 
the  streamwise  sectional  twisting  as  the  wing  statically  bends 
up.  Figure  2.1. a  shows  a  parallel  displacement  of  wing 
section.  Figure  2.1.b  has  an  induced  negative  angle  of  attack 
which  tends  to  restore  the  wing  deformation  in  the  air  flow  to 
the  original  position.  Figure  2.1.c  shows  the  induced  positive 
angle  of  attack  and  its  divergence  tendancy  with  airflow.  As  a 
result,  the  ply  orientation  of  Figure  2.1.b  will  undo  the 
aeroelastic  deformation  and  consequently  eliminate  or  delay  the 
aeroelastic  instability.  It  should  be  noted  that  the  dynamic 
type  behavior,  such  as  wing  flutter,  is  much  more  complicated 
than  the  simple  static  mechanism  illustrated  above.  This 
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simple  static  example,  however,  should  provide  understanding  of 
the  underlying  mechanism  and  initial  trends  of  the  aeroelastic 
tailoring . 


2.2  Aeroelastic  Formulation  and  Analysis 


The  flutter  analysis  of  an  aircraft  in  steady  flight  is 
formulated  by  the  Rayleigh-Ritz  approximation.  The  development 
of  the  mathematical  model  governing  the  uncoupled  longitudinal 
motions  considers  small  disturbances  from  a  steady  rectilinear 
flight.  The  small  disturbances  include  the  vehicle  rigid  body 
motions  and  elastic  deformations,  and  are  described  in  an  axis 
system  moving  steadily  with  the  vehicle  in  still  air.  For 
steady  level  flight,  the  axis  system  can  be  considered  fixed  in 
space  while  the  air  is  flowing  in  the  opposite  direction  at 
constant  speed,  and  is  identical  to  the  tunnel  support  axis 
system  as  shown  in  Figure  2.2.  To  isolate  and  emphasize  the 
effect  of  wing  aeroelastic  behavior,  the  analysis  hypothesizes 
a  rigid  aircraft  with  the  only  flexible  component  being  the 
wing.  The  governing  equation  of  the  aeroelastic  system  is 
derived  by  the  energy  method,  i.e.,  a  Rayleigh-Ritz  type 
formulation.  This  method  has  the  advantage  of  using  an  axis 
system  which  does  not  coincide  with  the  center  of  mass  of  the 
deformed  configuration. 


a)  Structural  Modeling 
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Any  vertical  displacement  of  the  vehicle  with  respect  to  a 
space-fixed  horizontal  plane  is  approximated  by  a  superposition 
of  finite  number  of  assumed  modes,  i.e., 


w  ( x , y , t )  -  j  Yjjx.y)  q- ( t ) 


(2.2) 


-  E  ♦i(*)*i(y)qi(t) 

i-1  1  1  1 


where  N  is  the  number  of  assumed  modes.  The  assumed  modes, 
ri(x,y),  include  aircraft  vertical  translation  (plunging)  and 
pitching  displacement  about  the  undeformed  aircraft  center  of 
gravity,  and  elastic  wing  component  modes.  The  elastic  wing 
component  modes,  conveniently  expressed  in  the  wing  coordinate 
system,  represent  wing  deformations  relative  to  the  wing  root 
position  for  the  rigid  fuselage  assumption.  The  degree  of 
accuracy  provided  by  the  Rayleigh-  Ritz  type  analysis  depends 
entirely  upon  the  choice  and  number  of  modes.  For  example,  the 
use  of  unrestrained  normal  modes  will  give  more  accurate 
results  and  fast  convergence.  It  is,  however,  structurally  and 
aerodynamically  convenient  to  use  the  uncoupled  beam-type 
bending  and  torsional  modes  for  the  equivalent  plate  model.  In 
addition.  Ref.  30  has  revealed  the  importance  of  the  additional 
chordwise  mode  in  modeling  vibration  analysis  of  an  anisotropic 
plate.  Using  the  same  assumed  wing  component  modes  in  Refs.  9 
and  10,  the  seven  assumed  modes  are 
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where  -  1.8751,  -  0.7341,  e2  -  4.6941  and  a2  -  1.0185. 

The  derivation  of  governing  differential  equations 
involves  the  evaluation  of  kinetic  energy,  T,  and  potential 
energy,  U,  and  the  use  of  Lagrange's  equation,  i.e., 


dr 
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where  L 

-  T  - 

U  and 

for  i-1  to  7 


(2.4) 


via  the  virtual  work  expression.  The  kinetic  energy  of  an 
unrestrained  aircraft  is 


T  -  J  II  "A  “2  dX  d* 


(2.5) 


where  (’)  -  d/dt  and  mA  is  the  mass  per  area  in  x-y  plane. 
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Because  of  the  zero  potential  energy  associated  with  the 
aircraft  rigid  body  motions,  the  total  potential  energy 
contains  only  the  strain  energy  of  the  elastic  wing.  In 
practice,  only  the  out  of  plane  deformation  is  relevant  to  the 
aerodynamic  operator.  Therefore,  only  the  bending  strain 
energy  will  be  evaluated  for  a  symmetric  laminate  plate  wing, 
i  .  e . , 


u  ■  i  [DUW-J5  ♦2D12"'i;w-yy+D22"'Y  ♦' «D16W' i*”-  5y 


wing 


(2.6) 
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where  w,-  denotes  3w/3x  and  so  on,  and  the  D. .'s  are  laminate 

X  1  ] 

bending  stiffness  terms  derived  from  classical  laminate  plate 
theory  (29].  Substituting  Eq.  (2.2)  into  (2.5)  and  (2.6),  one 
obtains  the  kinetic  and  potential  energy  in  terms  of 
generalized  coordinates  as 
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Appendix  A  details  the  inertia  and  stiffness  coefficients.  In 
analyzing  the  tunnel  restrained  configuration,  the  additional 
potential  energy  terms  are  added  as 


U  -  I  W  *  I  k2*22 


(2.9) 


where  and  K2  are  the  support  plunge  and  pitch  spring 
constants.  Placing  Eq.  (2.7)  and  (2.9)  into  Eq.  (2.4)  results 
in  a  set  of  differential  equations  of  motion  of  the  flexible 
aircraft  as 


M  {q}  +  K  {q}  -  {Q} 


(2.10) 


where  { q }  -  (q,,  q, ,  q,f  q4,  qg,  qg ,  q? ) 


(b)  Transient  Aerodynamic  Modeling 


For  the  purposes  of  the  present  study,  it  is  intended  to 
develop  an  accurate  and  relatively  easy  flutter  analysis  tool 
to  interpret  and  correlate  the  experimental  aeroelastic 
behavior  of  a  three-dimensional  vehicle  configuration.  The 
two-dimensional  incompressible  aerodynamic  theory,  including 
camber  oscillation  [31],  is  hence  chosen  even  though  there 
exists  advanced  but  computationally  complicated  unsteady 
aerodynamic  theory.  To  compensate  the  shortcoming  of  using 
two-dimensional  aerodynamic  tool  for  a  complicated  vehicle 
configuration,  two  modifications  are  made  to  the  proposed 
aerodynamic  tool.  One  is  the  so-called  "modified  strip  theory” 
[32]  which  accounts  for  finite  span  effect  by  using  the 
spanwise  lift-curve-slope  distribution.  The  static 
lift-curve-slope  distribution  is  obtained  by  Weissinger 
L-method  with  constant  angle  of  attack  across  the  span  [33]. 
This  approximate  procedure  uses  a  lifting-line  type  of 
representation  of  vortex  sheet,  but  satisfies  downwash  boundary 
condition  at  the  three-quarter  chord  point  of  each  spanwise 
station.  This  implies  a  lift-curve-slope  of  2n  in  the 
specializing  case  of  a  two-dimensional  airfoil.  In  Ref.  33, 
charts  of  the  necessary  influence  coef f icients ,  which 
facilitate  wing  loading  calculation,  are  presented.  Althugh 
the  procedure  in  Ref.  33  involves  the  evaluation  of  wing 
loading  at  only  seven  spanwise  stations,  Multhopp  integration 
points,  the  resulting  accuracy  is  considered  adequate  for  the 


present  purpose,  and  the  method  is  used  because  of  its 
simplicity. 

For  the  present  aircraft  model,  30°  forward  swept  wing 
with  aspect  ratio  of  seven,  the  symmetric  lift-curve-slope  are 
given  at  four  Multhopp  integration  station  (semi-span)  as 

Cn  -  2.20  at  92%  semi-span 
'a 

CQ  -  3.63  at  70%  semi-span 
*•  a 

C n  «  4.54  at  38%  semi-span 
a 

C n  -  5.74  at  fuselage  centerline 

*  a 

The  spanwise  lift-curve-slope  distribution  used  in  the  flutter 
analysis  is  then  interpolated  at  ten  equally  divided  wing  spans 
as  shown  in  Figure  2.3  which  gives  an  average  value  of  4.07 
across  the  full-span. 

The  spanwise  distributions  of  the  lift  and  moment  on  the 
deformed  wing  are  then  found  by  using  the  aforementioned  values 
of  static  sectional  lift-curve-slope  in  conjunction  with  the 
effective  angle  of  attack  distribution  resulting  from  wing 
oscillation.  As  suggested  in  Ref.  32,  the  proposed 
modification  should  be  applied  only  to  cases  involving  low  to 
moderate  reduced  frequencies.  The  examination  of  the  results 
of  flutter  analysis  is  neccessary  to  validate  the  assumption. 

The  other  modification  is  intended  to  account  for  the 
canard/wing/fuselage  interaction  by  using  the  "rigid"  aircraft 
lift  and  moment  characteristics  measured  experimentally.  The 
total  lift  and  moment  associated  with  the  rigid  body  motions 
are  then  evaluated  by  using  measured  aerodynamic  derivatives, 
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CL  ,  CM  ,  and  CM#/  as 
a  <*  a 


LA-  9SwCL  («2-«l/V)  “  ^SwcCM  VV 


(2.11 


Ma"  «SwcCh  ^2"<5l/V)+  <JSwc  CM-VV 


where  q  is  the  dynamic  pressure,  Sw  and  c  are  the  total  wing 
area  and  mean  aerodynamic  chord  length  respectively. 

The  transient  aerodynamic  coefficients  are  then  derived 
from  the  expression  of  virtual  work  due  to  pressure  disturbance 
as 


SWg  «  JJ  Ap  Sw(x,y)  dx  dy 
A/C 


(2.12) 


In  light  of  the  proposed  semi-analytical  aerodynamic  tool,  it 
is  convenient  to  decompose  the  virtual  work  contribution  into 
"body"  and  "wing"  parts  as 


SW. 


-  Ap  Sw  dxdy 
A/C 


■  Lbody  5ql  *  "body  Sq2  *  J  4P  5”  dxdlr 


(2.13) 


wing 


where  "body"  denotes  contributions  from  canard/fuselage 
contribution.  Incorporating  wing  airloads  described  by  the 
modified  streamwise-str ip  theory,  the  virtual  work  expression 
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is  then 


'body 


6ql  +  Mbody  Sq2  +  J  (Lih  +  M5ot  +  N5*>  dx 


wing 


(2.14) 


- 


where  L,  M  and  N  are  sectional  forces  and  moments,  and  h,  a, 
and  l  represent  sectional  deformations  defined  at  sectional 
mid-chord  position.  For  a  Rayleigh-Ritz  type  analysis,  the 
spanwise  variation  of  sectional  deformations  are  defined  and 
expressed  in  terms  of  the  assumed  modes  as, 


ql  +  (yw“xsinA)q2  +  *3q3  +  *5q5  "  E  niqi 

i  *  1 


-  q2  -  +  ^5)sinA  +  ♦4c34  +  ^6c36)cosA 


(2.15) 


i-l  1  1 

*7q7  *  E^i^i 


Since  the  streaowise  section  is  not  normal  to  the  structural 
span  of  a  swept  wing,  the  sectional  twist  and  camber  change  are 
closely  approximated  by  the  assumed  modes  in  Eq.  (2.15).  In 
here,  A  and  yw  are  both  negative  values  for  the  forward  swept 
wing  aircraft  shown  in  Figure  2.2. 

In  studying  the  interaction  between  elastic  and  rigid  body 
motions,  not  only  the  instability  onset  but  also  the  aircraft 
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flying  quality  are  of  great  concern  [34J.  The  aircraft  flying 
quality  is  characterized  by  its  transient  behavior  such  as 
decay  rate,  natural  modes  and  frequencies.  To  describe  these 
transient  characteristics,  airloads  due  to  arbitrary  wing 
motions,  instead  of  harmonic  airloads,  must  be  used.  One  such 
method,  in  two-dimensional  incompressible  flow,  is  to 
approximate  the  generalized  Theodorsen  function,  C(p),  in  the 
Laplace  domain  by  a  rational  polynominal, 


C(p) 


biP  +  b0 
P  +  b0 


(2.16) 


where  bj^  «  0.55,  bg  -  0.15,  p  is  the  dimensionless  Laplace 
variable  p  •  pb/V  ,  and  p  -  a  +  jw.  The  above  single  pole 
approximation  gives  a  reasonable  fit  to  the  Theodorsen  function 
over  the  reduced  frequency  k  -  bw/V  from  0  to  1  as  shown  in 
Figure  2.4.  It  should  be  noted  that  the  high  frequency 
asymptote  of  the  above  approximation  approaches  0.55  instead  of 
the  exact  value  of  0.50  which  suggests  some  inaccuracy  of  the 
current  approximation  outside  the  low  reduced  frequency  range. 
For  the  present  study,  the  reduced  frequency  of  body  freedom 
flutter  ranges  from  0.02  to  0.05.  The  transformed  airloads,  as 
detailed  in  Appendix  B,  can  then  be  described  in  the  Laplace 
domain  as 
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Up)  -  Cg  pv  b 

*•  a 

{[B2AP2  +  B1AP  +  B0A  +  B3A  P/(P+bO)]  h/b 

+[B2BP2  +  B1BP  +  B0B  *  B3B  P/(P+bO)]  a 

+[B2CP2  +  B1CP  +  B0C  +  B3C  P/(P+bO)]  */b 

M(p)  -  Cg  pV^b^  |[*‘]  +  [•*]  <*  +  [  ’  *  j  £/b  | 

N ( p )  -  C £  pV^b  |[‘‘]  h/b  +  [•*]  <*  +  [**]  t/b  | 


where  the  B2a'  b1a'  B0a  anc*  B3A  etc  rePresent  th®  aerodynamic 

virtual  mass,  damping,  stiffness  and  lag  terms  respectively, 

and  Cg  is  the  streamwise  sectional  lift-curve-slope.  Placing 
a 

Eq.  (2.17)  into  (2.14),  introducing  h,  a  and  i  from  Eq.  (2.15), 
and  defining  new  augmented  state  variables  yi  such  that 


y  (p)  -  — -  q Up) 

P  +  b0 


for  i  -  1  to  7 


(2.18 


allows  one  to  take  the  inverse  Laplace  transform  and  write  the 
generalized  aerodynamic  forces  in  transient  form  as 


■  Z  IqA-  q.  +  QB  •  q-  +  QC.  q.  +  Qb  .  y  . ) 
j_l v  i 3  43  43  13  43  13  y3^ 


(2.19) 
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*i  -  q: 


for  i  -  1  to  7 


(2.20) 
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where  Q^'s  are  the  aerodynamic  coefficients  detailed  in 
Appendix  B. 
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( c )  Aeroelastic  Analysis 

The  augmented  states,  introduced  in  the  transient  airloads 
modeling,  increase  the  order  of  the  system  but  allow  the 
convenience  of  linear  eigenvalue  analysis  of  the  constant 
coefficient  ordinary  differential  equations.  Placing 
Eq.  (2.19)  into  (2.10)  and  incorporating  the  additional 
governing  equations  of  aerodynamic  lags,  the  resulting 
aeroelastic  equation  of  motion  are  then  expressed  in  the  state 
variable  form  as 


10  0 

o  (m-qa)  o 
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(qc-k)  (Qb-c)  qd 

/  \ 
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0  I  -H 
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\ 

(2.21) 


where  I  is  the  identity  matrix  and  H  ■  bgV/b  I.  This  yields  a 
set  of  twenty-one  first  order  differential  equations.  For 
given  values  of  dynamic  pressure,  the  transient  behavior  and 
stability  characteristics  are  determined  by  setting  q  -  qe^,  y 
-  ye^>t,  and  finding  the  eigenvalues  p  of  Eq.  (2.21).  The 
flutter  and  divergence  boundaries  are  determined  by  solving  the 
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successive  eigenvalue  problems  until  the  instabilities  occur. 
The  transient  aeroelastic  behavior  is  represented  by  the 
successive  eigenvalues  and  eigenvectors.  If  the  assumed 
elastic  modes  are  set  equal  to  zero,  the  above  stability 
analysis  reduces  identically  to  the  rigid  aircraft  dynamic 
stability  analysis  in  the  space-fixed  axis  system. 

Also,  it  is  to  be  noted  that  since  modern  control  theory 
is  well  developed  for  systems  described  by  finite  order, 
constant  coefficient  ordinary  differential  equations,  the  other 
motivation  for  using  transient  airload  expressions  is  the 
potential  of  active  control  design  interaction  with  the 
aeroelstic  analysis. 

For  comparison  purposes,  the  cantilever  wing  flutter  and 
divergence  analyses  can  be  performed  by  setting  rigid  body 
modes  to  zero.  For  wing  divergence  analysis,  the  finite  span 
effect  will  be  accounted  for  by  applying  an  overall  correction 
factor  suggested  in  Ref.  35  as 


d_  c  -  a  AR 

do  L  o  AR  +  4 


(2.22) 


where  8q  is  the  two-dimensional  lift-curve-slope  for 
incompressible  flow,  2n,  and  AR  is  the  aspect  ratio,  2-£/c  of 
the  wing.  This  is  equivalent  to  modifying  the  wing  divergence 
speed  VD,  found  from  the  2n  lift-curve-slope,  by  the  square 
root  of  the  overall  correction  factor  as 
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AR  +  4 
AR 


<V2n 


(2.23) 


For  the  cantilever  wing  flutter  analysis,  no  finite  span 
correction  is  used.  This  is  because  of  the  close  resemblance 
to  two-dimensional  airload  distribution  due  to  the  high 
frequency  wing  oscillatory  motion  at  flutter  condition. 

2 . 3  Wind  Tunnel  Support  Stability  Analysis 

In  a  "free  flight"  wind  tunnel  test,  the  support  system 
not  only  physically  restrains  the  model  freedoms,  but  may  also 
influence  the  overall  dynamic  characteristics  and  corresponding 
testing  results.  Without  careful  consideration  of  the  support 
system  dynamics,  the  testing  results  may  be  contaminated  and  an 
additional  support  instability  may  exist.  For  such  "free 
flight"  wind  tunnel  tests,  various  remedies  have  been  developed 
as  "ad  hoc"  compensation  schemes  with  some  degree  of  success, 
but  no  clear  analytical  explanation  has  been  given. 

A  simplified  longitudinal  stability  analysis  model  is 
shown  in  Figure  2.5.  The  aircraft  model  is  considered  to  be 
"rigid"  in  this  analysis.  This  support  configuration  provides 
freedoms  in  vertical  translation,  pitch  and  yaw  motions.  with 
the  assumptions  of  small  disturbances,  the  governing  equations 
for  the  longitudinal  mode  are 
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(2.24) 


where  M  and  I  are  model  mass  and  pitch  moment  of  inertia,  and 
,  K2  are  the  support  plunging  and  pitching  stiffness 
constants.  The  yaw  motion  is  uncoupled  here  from  the 
longitudinal  motion,  i.e.,  plunge  and  pitch.  Placing 
Eq.  (2.11),  La  and  MA,  into  (2.24),  the  equation  of  motion 
becomes 
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where  (°)  ■  d/dr  and  wr  is  the  reference  frequency.  The  modal 
viscous  damping  ratios,  and  C2,  are  introduced  into 
Eq.  (2.27)  to  account  for  the  dissipating  force  and  moment  in 
the  support  system.  Rearrange  Eq.  (2.27)  into  matrix  form  as 
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where 
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C12  “  JJJ  CM  9 
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For  a  linear  stability  problem,  assuming  solution  in  a  form  of 
e^T  yields  the  stability  characteristic  equation  as 


(p  +C^P+K^  )  (  cjt2P+Ki2  ^ 

(C21P+K21)  (P2+C22P+K22) 


(2.30) 


E4p  +  E3p  +  E2p  +  Ej_p  +  EQ  -  0 


(2.31) 
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C11  +  C22 

K11  +  K22  +  C11C22  '  C12C21  (2.32) 

K11C22  +  K22C11  '  K12C21 

K11K22 


The  support  system  stability  can  then  be  determined  by 

examining  the  real  part  of  all  the  roots  of  Cq.  (2.31).  For  a 
fourth  order  polynominal  with  multiple  parameters  embedded  in 
its  coefficients,  one  can  determine  the  roots  directly  for 
given  values  of  parameters.  By  doing  so,  however,  the 

parameter  dependency  and  the  intrinsic  mechanism  of  stability 
may  be  lost. 

To  obtain  a  closed  form  stability  criterion  without  losing 

its  generality,  the  equation  of  motion  is  simplified  by 

ignoring  the  viscous  damping  ratios,  and  £2*  and  aerodynamic 

derivative  CM> .  By  determining  the  sign  of  the  real  part  of 
a 

all  the  roots,  the  Routh-Hurwitz  criterion  is  useful  in  giving 

a  stability  condition  in  terms  of  polynominal  coefficients  and 

system  parameters.  The  real  part  of  all  the  roots  are 

negative,  hence  asymptotically  stable,  if  and  only  if  all  the 

coefficients  E^,  E^,  E2,  E^,  Eg  and  the  discriminant 

2  2 

D  -  Ej^Ej  -  EqEj  -  E^E^  are  positive.  Placing  Eq.  (2.29) 
into  (2.32),  the  coefficients  E.  and  discriminant  D  become 
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For  a  statically  stable  aircraft,  CL  is  positive  while  CM  is 

a  a 

negative.  These  coefficients  are  then  all  positive,  and  the 
condition  of  positiveness  for  the  discriminant  D  provides  a 


constraint  equation  for  the  support  stiffness  required  for 


stability,  namely 


(2.34) 


Except  for  a  very  lightweight  model,  the  second  term, 


CM  /In,  of  Eq.  (2.34)  is  usually  small  and  can  be  neglected. 


Therefore,  Eq.  (2.34)  simply  states  that  the  support  plunging 


frequency  must  be  lower  than  the  pitching  frequency  in  order  to 


avoid  the  unwanted  support  instability.  This  is  an  important 


design  criterion  for  the  support  stiffness  of  a  "free  flight" 
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wind  tunnel 

test 

setup. 

It 

should  be 

noted  that  the  above 

analysis  is 

not 

an  "ad 

hoc" 

solution 

to  the  vertical  rod 

support  system  used  herein.  The  support  stiffness  can  arise 
from  any  non-mechanical  means  such  as  cable  restraining  forces 
of  a  two  cable  mount  system.  The  above  criterion  of  support 
stiffness  should  be  applicable  to  other  type  of  wind  tunnel 
support  system. 

By  excluding  ,  C 2  and  CM>,  Eq.  (2.34)  gives  a  sufficient 

a 

stability  criterion  which  is  independent  of  the  airspeed. 

Since  ,  C2  and  CM<  are  all  positive  damping  terms  in  diagonal 

a 

positions,  and  C22'  t*ie  inclusion  of  these  damping  terms 

will  generally  stabilize  the  system.  As  a  result,  there  may 
exist  a  stability  boundary  or  region  as  a  function  of  airspeed 
in  addition  to  the  stability  criterion  of  Eq.  ^2.34) .  The 
complete  stability  solution  can  only  be  determined  by  the 
successive  polynominal  solutions  of  Eq.  (2.31).  Figure  2.6 
demonstrates  the  effect  of  plunging  damping  ratio  in  a  complete 
stability  analysis.  With  a  finite  amount  of  plunging  damping 
ratio,  the  support  system  will  be  eventually  stabilized. 

The  physical  interpretation  of  support  instability  arises 
from  the  interaction  between  the  plunge  and  pitch  modes  of  a 
"rigid"  model.  Since  the  pitch  frequency  increases  with 
airspeed  due  to  aerodynamic  stiffening  while  the  plunge 
frequency  is  independent  of  the  airspeed,  there  will  be  a 
frequency  coalescence  and  mode  transfer  if  the  pitch  frequency 
in  still  air  is  below  the  plunge  frequency.  Figure  2.7 
illustrates  this  interaction  phenomenon  on  an  w-V  plot. 
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Fig.  2.6  Plunge  Damping  Effect  on  Support  Instability 
(Root  Locus  Plot,  Rigid  Aircraft  Model) 
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For  the  purposes  of  this  investigation,  a  generic, 
full-span,  forward  swept  wing  aircraft  model  was  constructed. 
Figure  3.1  shows  the  layout  of  the  aircraft  model  and 
Figure  3.2  shows  picture  of  the  installed  model  on  wind  tunnel 
site.  The  principal  components  are  (1)  the  inte rchangable , 
30-degree  forward  swept,  graphite/epoxy  flat  plate  wings,  (2)  a 
sturdy  fuselage  of  balsa  wood,  and  (3)  a  rigid,  all  moving 
canard  surface. 

The  wing  planform  was  chosen  to  be  the  same  as  those  in 
Refs.  9  and  10  such  that  the  present  study  complemented  the 
previous  cantilever  wing  studies  and  isolated  the  effect  of 
rigid  body  freedoms.  Instead  of  using  semi-span  wing  models  as 
in  Refs.  9  and  10,  full-span  wing  models  were  fabricated  by 
bonding  two  305  x  76  ram  (12  x  3  inch)  laminate  plate  onto  the 
fiberglass  center  plate  in  a  30°  forward  swept  configuration  as 
shown  in  Figure  3.3.  The  laminate  plate  wings  were  made  up  of 
Hercules  AS4/3501-6  graphite/epoxy  which  is  supplied  in  305  mm 
wide,  resin  impregnated  tape  ("prepreg").  Material  properties 


for  this  material  are  indicated  in  Table  3.1. 

The  manufacturing  of  laminates  followed  the  standard 


Table  3.1  Graphite/Epoxy  Material  Properties 


Hercules  AS4/3501-6 


Property 

Extensional 

Flexural 

El  (Gpa) 

143. 

106  . 

Et  (Gpa) 

9.7 

7.9 

\)  . 

.  3 

.  3 

LT 

Glt  (Gpa) 

4.9 

4.9 

Measured  Ply  Thickness  -  0.131  mm 


1  obtained  from  vibration  tests 

2  AS1/3501-6  property 
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procedures  developed  at  TELAC  over  the  years  (36)  except  for 
the  placement  of  paper  bleeder  for  curing.  This  difference 
will  be  detailed  in  the  following  description  of  laminate 
manufacturing.  The  prepreg  was  first  removed  from  its  storage 
at  -18°C  and  cut  to  size  using  Stanley  knives.  Templates  were 
used  to  cut  all  the  desired  angle  plies.  The  plies  were  then 
laid-up  in  350  mm  by  305  mm  laminates  and  covered  with  peel-ply 
on  both  sides.  The  laminates  were  then  placed  on  an  aluminum 
cure  plate  sprayed  with  mold  release  and  covered  with 
non-porous  teflon.  Instead  of  placing  three  sheets  of  paper 
bleeder  only  on  top  of  each  6-ply  laminate,  two  sheets  of  paper 
bleeder  were  placed  on  each  side  of  all  the  laminates  for 
curing.  Due  to  the  manufacturing  residual  strain,  previous 
symmetric,  thin  laminate  often  exhibited  some  degree  of  warping 
and  unintended  spanwise  twisting.  The  purpose  of  this 
modification  in  manufacturing  procedure  was  to  minimize  or 
eliminate  any  symmetric  laminate  warping  by  using  a  "symmetric" 
curing  setup.  Other  necessary  curing  materials  were  stacked  as 
well  and  the  complete  setup  is  illustrated  in  Figure  3.4. 

The  cure  cycle  used  for  the  3501-6  epoxy  was  a  two-step 
process  composed  of  a  one-hour  flow  stage  at  116°C  and  two-hour 
set  stage  at  176°C.  Full  vacuum  was  applied  to  the  laminates 
and  a  0.59  Mpa  (85  psi)  pressure  was  maintained  inside  the 
autoclave  during  the  entire  cure  (see  Figure  3.5).  After  that, 
the  laminates  were  postcured  for  eight  hours  at  176°C  in  an 
oven  with  no  pressure  or  vacuum. 

Using  a  milling  machine  equipped  with  a  diamond  grit  blade 
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Fig.  3.5  AS4/3  501-6  Graphite/Epoxy  Cure  Cycle 


<r  «-* 


and  water  cooling,  three  wing  strips  (350  x  76.2  mm)  were  cut 
out  of  each  laminate.  From  the  same  laminate  two  wing  strips 
were  selected  to  form  the  full-span  wing  by  putting  one  of  the 
wing  strips  upside  down  such  that  the  ply  orientation  would  be 
symmetric  about  the  fuselage  center  line.  The  wing  strips  were 
then  sandwiched  bonded  to  between  two  Scotchply  1003 
fiberglass,  V-shaped  center  plates.  The  epoxy  used  for  bonding 
was  an  Epon  815  resin  and  V40  hardener  supplied  by  the 
Mi 1 le r-Stephenson  Company.  The  mixing  ratio  used  was  60:40, 
resin  to  hardener  by  weight.  The  cure  used  was  room 
temperature  cure  last  for  seven  days.  During  the  cure,  the 
wing/center  plate  assemblies  were  firmly  held  together  by  five 
C-clamps . 

The  effects  of  aeroelastic  tailoring  were  demonstrated  by 
a  set  of  four  wings  with  bending-twi sting  coupling  in  both 
favorable  and  adverse  fashions,  namely,  [02/90  1,  *  ( 1 52/0 ] s , 
[302/0]g,  and  [-152/0]g  layups.  The  use  of  symmetric  but 
unbalanced  laminates  provided  a  possible  wide  range  of 
bending-twisting  coupling.  The  bending-twi sting  coupling 
stiffness  ratios,  D^g//  D^  Dgg  and  the  bending/twisting 
stiffness  ratios,  Dn/°5g»  of  all  the  specimens  are  tabulated 
in  Table  3.2.  For  all  wings,  the  zero  degree  ply  was  used  as 
the  inner  ply  to  strengthen  the  intended  "unidirectional" 
laminate  except  for  the  (02/90)s  wing  in  which  the  inner  ply 
angle  was  90°  to  avoid  true  unidirectional  laminate. 

The  fuselage  section  was  a  sturdy  box  structure  design 
with  9.5  mm  (3/8  inch)  thick  balsa  wood  and  a  pine  wood  nose 
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3.2  Coupling  Stiffness  and  Bending/Torsion  Stiffness 
Ratios  of  Composite  Wings 


D16//Dll  °66 

Dll/D66 
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cone.  The  wind  tunnel  support  ran  through  the  open  holes  of 
the  fuselage,  through  the  support  interface,  i.e., 
pitching/plunging  bearing  assembly,  inside  the  fuselage.  Two 
access  doors  were  designed  on  top  of  the  fuselage  for  easy 
assembly.  The  wing  mount  was  underneath  the  rear  fuselage 
section.  Five  built-in  T-nuts  were  used  to  assure  the  mount 
rigidity.  The  wing  mount  section  had  a  step  recess  so  that  the 
externally  mounted  wing  would  be  flush  with  the  fuselage  bottom 
surface . 

The  airfoil  shaped  canard  surfaces  were  mounted  on  a 
single  shaft  running  through  the  mid-height  of  the  fuselage 
sidewalls.  Two  teflon  collars  were  used  for  the  shaft/fuselage 
intersections.  The  canard  was  positioned  by  two  clearance 
wedges . 

The  wing  gage  wires  were  directed  inside  of  the  fuselage 
through  a  passage  hole,  and  lead  outside  of  the  fuselage  along 
the  suspension  spring.  All  the  wires  were  slackly  attached  to 
the  suspension  spring  to  minimize  the  wire/model  interference. 
It  was  particularly  advantageous  to  handle  the  wires  near  the 
pitching  axle  so  that  the  pitching  stiffness  was  less 
interferred  with. 

{ b )  Wind  tunnel  support  system 

Wind  tunnel  investigation  of  the  body  freedom  flutter 
requires  that  the  dynamic  behavior  of  a  freely  flying  model  be 
simulated.  In  the  wind  tunnel  study,  however,  the  model 
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generally  must  be  somehow  restrained  and  the  behavior  may  be 
altered  to  a  certain  degree  by  the  support  system-  without 
careful  consideration  of  the  dynamics  of  the  overall  system,  an 
inherent  support  instability  may  exist.  Various  model  support 
systems  [37]  were  studied  for  the  MIT  low  speed  acoustic  wind 
tunnel  which  has  a  1.5  x  2.3  m  (5  x  7.5  ft)  free  jet  test 
section . 

The  general  requirements  were  to  have  a  soft,  stable  and 
simple  model  support  system  with  negligible  moving  mass  and  low 
aerodynamic  interference.  Because  of  the  physical  limitation 
of  minimum  model  weight,  the  preselected  wing  planform  was  not 
able  to  fly  the  model  completely.  This  overweight  factor, 
together  with  the  requirement  of  soft  support,  lead  to  the 
selection  of  a  vertical  rod  support  system.  The  vertical  rod 
support  system  was  originally  developed  by  Boeing  Airplane 
Company  and  has  been  successfully  used  in  low  speed  wind  tunnel 
tests  [38].  The  modified  support  system  for  the  MIT  acoustic 
wind  tunnel  included 

(1)  a  1.83  m  (6  ft)  long,  vertically  installed,  Case 
hardened  steel  rod  with  a  19  mm  (3/4  inch)  diameter, 

(2)  a  Thomsoa  super-12  linear  bearing  sliding  on  the  rod 
and  attached  to  the  fuselage  through  a  pitch  bearing  mount 
located  at  the  model  center  of  gravity, 

(3)  a  soft  helical  spring  attached  to  the  linear  bearing 
which  suspended  the  model  from  the  top  of  the  vertical  rod, 

(4)  a  pair  of  linear  springs  attached  to  the  pitch  gimbals 
at  forward  and  aft  extended  arms  which,  together  with  the 
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existing  model  pendulum  stiffness  provided  the  support  pitching 
stiffness,  and 

(5)  snubber  cables  attached  to  the  fuselage  at  the  nose 
and  tail  locations  which  lead  outside  the  test  section  to 
control  sticks. 

With  this  test  setup,  the  model  was  softly  restrained  in 
the  plunge,  pitch  and  yaw  motions,  but  totally  restrained  in 
roll,  lateral  and  fore-and-aft  motions  except  for  slight  free 
play  between  the  linear  bearing  and  rod.  Figures  3.6  and  3.7 
show  the  wind  tunnel  setup.  Details  of  the  model  properties 
and  model  support  dynamic  characteristics  are  given  in 
Table  3.3. 

Since  the  linear  bearing  could  develop  significant  amount 
of  friction  from  misalignment  or  applied  moment,  the  properly 
assembled  linear  bearing  and  moment  free  condition  became 
extremely  important.  In  the  model/support  assembly,  however, 
the  helical  suspension  spring  always  applied  a  moment  on  the 
bearing  due  to  the  eccentric  force  at  the  spring  periphery  The 
spring/model  attachment  was  therefore  arranged  to  cause  only 
model  roll  unbalance  but  no  pitch  unbalance.  This  roll 
unbalance  was  subsequently  counterbalanced  by  a  roll  balanced 
weight.  The  so-arranged  model/support  assembly  was  well 
balanced  in  pitch  and  coll  before  putting  in  the  linear 
bearing . 


MIT  ACOUSTIC  WIND  TUNNEL 
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Fig.  3.6  Wind  Tunnel  Setup 
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Table  3.3  Aircraft  Model  and 

Support  Properties 

-3 

V 

Mass 

0.956  Kg 

«  V 

2 

S' 

Pitch  Inertia 

0.0178  Kg-m 

Wing  Mass 

0.054  Kg 

Wing  Span 

0.612  m 

Wing  Swept  Angle 

-30° 

“j 

Wing  Area 

0.0538  m2 

■-'J 

>  > 

Wing  Aspect  Ratio 

7 

;■ 

Canard  Area 

0.0126  m2 

fm 

Fuselage  Length 

0.8  m 

f'm 

Model  Plunge  Freq. 

0.63  Hz 

Model  Pitch  Freq. 

(original ) 

0.20  Hz 

•  m 

«  , 

Model  Pitch  Freq. 

( modified ) 

0.85  Hz 
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(c)  Instrumentation 


The  model  was  instrumented  to  monitor  the  aeroelastic 
behavior  including  the  model  rigid  body  motions  and  wing 
elastic  deformation.  The  wing  root  bending  and  torsional 
strain  gages,  on  both  sides  of  wings,  monitored  the  wing 
motions.  Figure  3.8  shows  the  arrangement  of  strain  gages. 
The  bending  gages,  M-M  EA-13-187BB-120 ,  were  mounted  on  both 
faces  of  the  wing  to  form  a  half  Wheatstone  bridge.  To  monitor 
the  torsion  response,  the  2-element  90°  rosette  gages,  M-M 
EA-06-250TK-120 ,  were  used,  and  a  full  Wheatstone  bridge  was 
formed  by  rosette  gages  on  both  faces  of  the  wing. 

A  low  friction,  conductive  film  type  potentiometer  was 
used  to  monitored  the  model's  pitch  motion.  This  potentiometer 
was  snugly  connected  to  the  pitching  axle  and  mounted  outside 
the  fuselage  as  shown  in  Figure  3.9. 

A  thin,  flexible  203  mm  (8  inch)  long,  0.51  mm  (.020  inch) 
thick,  cantilever  steel  beam  with  a  root  strain  gage, 
Figure  3.10,  was  mounted  near  the  top  end  of  the  vertical  rod 
to  record  model's  plunge  motion.  The  free  end  of  the 
cantilever  steel  beam  was  tied  to  near  the  second  coil  of  the 
suspension  spring.  The  movement  of  the  free  end  would  be 
fractions  of  the  model  plunge  height  as  long  as  the  spring 
deflection  was  within  the  linear  range. 

The  wing  bending  and  torsional  gage  signals  were  monitored 
on  oscilloscopes  and  recorded  on  a  strip  chart  recorder 
together  with  the  pitch  and  plunge  signals.  For  the  four 
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channel  recorder  used  in  the  present  testing,  only  a  selected 
side  of  the  wing  root  bending  and  torsional  gage  signals  plus 
the  pitch  and  plunge  signals  were  recorded.  The  left  side 
wing,  facing  the  tunnel  control  station,  was  usually  recorded 
while  the  right  side  wing  was  occasionally  recorded  for 
comparison  purposes.  During  the  tests,  video  movies  of  the 
model  responses  were  also  taken  for  future  replay  and 
demonstration . 

3 . 2  Vibration  Tests 

The  objectives  of  the  vibration  test  were  to  determine  the 
vibration  frequencies  of  all  wings,  and  to  verify  the  closely 
tuned  frequencies  of  the  both  sides  of  the  full-span  wing. 
Since  the  structure  geometry  was  simple  and  had  well  separated 
vibration  frequencies,  the  sine-dwell  testing  [39]  proved  to  be 
very  powerful  and  efficient.  The  wing  model  was  horizontally 
mounted  atop  a  Ling  420-1  shaker  which  was  driven  by  the 
amplified  sinusoidal  signal.  The  vibration  frequencies  were 
tuned  in  by  progressively  varying  excitation  frequency  and 
observing  the  steady  state  resonant  responses.  The  excitation 
and  response  signals  were  also  fed  into  an  oscilloscope  to 
observe  the  90°  phase  shift  during  resonance.  The  excitation 
power  was  properly  adjusted  to  observe  the  vibration  modes 
without  incurring  large  amplitude  vibration.  This  procedure 
was  repeated  for  both  sides  of  the  wing,  which  usually  did  not 
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resonate  at  exactly  the  same  frequency.  The  aforementioned 
vibration  tests  were  performed  up  to  the  third  vibration  mode 
for  all  wings. 


3.3  Aerodynamic  Tests 


The  objective  of  the  aerodynamic  tests  was  two-fold.  One 
was  to  determine  the  model  trimmed  flight  performance  over  the 
range  of  angle  of  attack.  The  other  was  to  measure  the 
longitudinal  aerodynamic  derivatives  of  the  rigid  model. 

The  model  center  of  gravity  was  designed  to  be  slightly 
below  the  pitching  axle  such  that  a  stable  pendulum  mode  was 
obtained.  To  maintain  a  common  reference  angle  within  a  test 
and  among  other  tests,  an  on-site  model  balancing  was  developed 
as  a  standard  pre-testing  procedure.  With  wind  off  and  model 
locked  in  plunge  motion,  the  model  was  balanced  in  the 
horizontal  position,  defined  by  a  level  and  tail  position  on  a 
retractable  scale,  by  placing  balance  weights  atop  the 
fuselage.  After  the  model  was  balanced,  the  pendulum  stiffness 
was  measured  by  placing  incremental  counterbalancing  weights,  1 
gm  each,  at  positions  six  inches  in  front  of  and  behind  the 
pitching  axle,  and  reading  the  corresponding  pitch  angles.  The 
pendulum  stiffness  was  then  determined  from  the  plot  of 
balancing  moment  versus  pitch  angle.  The  model  was  then 
unlocked  in  plunge  motion,  to  determine  the  plunge  spring 
characteristics . 
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Calibrations  of  the  model  pitch  angle  and  plunge  height 
were  performed  before  each  group  of  aerodynamic  tests.  First, 
friction  test  was  carried  out  by  examining  the  model  sticktion 
("dead  zone")  in  both  plunge  and  pitch  motions.  Within  this 
dead  zone,  the  restoring  force  or  moment  could  not  overcome  the 
friction  force  or  moment  such  that  the  model  could  stay 
wherever  it  was.  A  3.2  mm  (.125  inch)  plunge  dead  zone  was 
usually  regarded  as  acceptable,  otherwise  the  linear  bearing 
needed  to  be  reassembled.  With  wind  off,  the  pitch  dead  zone 
was  about  4°  to  5°  for  the  original  soft  pitch  support,  but  not 
noticeable  for  the  later  stiffer  support.  Since  the  hysteresis 
arose  mainly  from  the  pitch  potentiometer,  no  improvement  could 
be  gained  on  the  existing  instrumentation  to  reduce  the  pitch 
dead  zone  in  still  air.  Careful  attention  to  hysteresis  was 
taken  during  the  static  tests.  It  should  be  noted  that  the 
hysteresis  was  reduced  due  to  the  higher  aerodynamic  pitching 
stiffness  while  wind  was  on.  Model  pitch  angle  was  then 
calibrated  by  the  differential  reading  of  the  tail  positions  on 
a  retractable  scale.  Model  plunge  height  was  calibrated  by  a 
small  scale,  76.2  ram  (3  inch)  maximum,  for  aerodynamic  testing 
and  by  a  large  scale,  304.8  mm  (12  inch)  and  up,  for  flutter 
testing. 

The  aerodynamic  characteristics  of  the  "rigid"  model  were 
obtained  on  the  actual  flutter  model  in  the  wind  tunnel,  by 
setting  the  airspeed  at  5  m/sec  to  avoid  aeroelastic  effect. 
For  the  [02/90)5  wing,  this  testing  speed  was  only  6%  of  the 
wing  divergence  dynamic  pressure.  The  procedure  of  placing 
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counterbalancing  weights  was  similar  to  that  of  the  pendulum 

stiffness  measurement  except  the  tunnel  was  now  running  at  a 

low  speed.  The  balancing  weight  was  increased  at  2  gm  each 

step  until  the  model  was  stalled.  The  1  gra  increment  was  then 

used  to  fine  tune  the  stall  regions.  At  each  incremental  step, 

the  pure  aerodynamic  moment  was  calculated  by  substracting  the 

corresponding  pendulum  moment  from  the  counterbalancing 

moment.  Also,  the  net  lift  was  calculated  by  subtracting  the 

counterbalancing  weights  from  the  change  of  the  suspension 

spring  force.  The  lift  and  moment  coefficients  were  then 

calculated  from  CL  -  Lift/qSw  and  CM  -  Moment/qSwc,  and  plotted 

a  a 

against  the  pitch  angle  readings.  This  procedure  was  repeated 
for  several  canard  settings,  typical  values  being  0®,  +2.5®, 

and  +5®. 


3.4  Cantilever  Wing  Tests 


The  objectives  of  the  cantilever  wing  testing  were  to 
verify  the  analytical  predictions  and  to  obtain  the  baseline 
results  for  body  freedom  flutter  comparison.  The  cantilever 
wing  configuration  was  set  up  by  blocking  the  model's  plunging 
freedom  with  top  and  bottom  stops,  and  varying  the  wing  root 
angle  of  attack  (model  pitch  angle)  with  the  help  of  the 
snubber  cables.  The  testing  procedure  explored  the  linear  and 
nonlinear  flutter  and  divergence  boundaries.  The  testing 
variables  were  wind  tunnel  speed  and  wing  root  angle  of  attack 
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from  0  to  16  degrees.  For  each  wing,  the  wind  tunnel  speed  was 
initially  set  to  be  low  enough  such  that  no  flutter  or 
divergence  was  present  up  to  large  wing  root  angle  of  attack, 
i.e.,  15°  or  16°.  The  airspeed  was  then  increased  in  1  m/sec 
steps.  At  each  incremental  speed  setting,  the  flutter  point 
was  determined  by  gradually  increasing  the  wing  root  angle  of 
attack  in  both  the  positive  and  negative  sense,  until  flutter 
was  encountered.  This  procedure  was  repeated  until  the  flutter 
or  divergence  condition  could  not  be  avoided  for  any  setting  of 
the  wing  root  angle  of  attack.  In  other  words,  the  positive 
sense  of  critical  angle  of  attack  coincided  with  the  negative 
one . 


Freedom  Flutter  Tests 


The  testing  objectives  were  to  explore  the  aeroelastic 
behavior  and  instability  boundaries  of  the  "free"  flying  model, 
and  to  verify  the  analytical  model.  The  tunnel  airspeed  was 
the  only  primary  variable,  while  the  canard  setting  was  used  as 
a  secondary  variable  for  flight  trimming.  For  a  given  test,  a 
canard  setting  was  chosen  to  make  the  trimmed  flight  pitch 
angle  slightly  positive.  From  the  previous  aerodynamic  tests, 
the  model  pitch  angle  was  chosen  at  about  2°,  since  models 
flown  at  this  angle  had  the  most  linear  moment  characteristics 
and  a  positive  lift  to  raise  the  model  into  the  middle  of  the 


test  section. 
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For  the  given  canard  setting,  normally  2.5°,  the  flutter 
point  was  determined  by  step  increasing  the  tunnel  airspeed 
until  the  body  freedom  flutter  or  wing  flutter  was 
encountered.  For  the  purposes  of  correlating  the  transient 
behavior,  subcritical  measurements,  such  as  frequency  and  decay 
rate  of  an  impulse  response,  were  made  at  each  incremental 
speed  step.  The  impulse  responses  of  pitch  and  plunge  freedoms 
were  generated  by  two  seperate  "kicks"  at  the  tail  and  near  the 
center  of  gravity  locations.  The  wing  elastic  responses  were 
also  generated  by  kicking  at  the  wing  tip  location.  Because  of 
the  low  airspeed  and  open  test  section,  the  needed  "kick"  was 
easily  accomplished  by  using  a  1.5  m  (5  ft)  long  wooden  stick 
reaching  the  model  from  downstream.  In  general  this  technique 
was  simple  and  effective,  but  it  became  more  difficult  and 
dangerous  to  the  model  as  the  tunnel  airspeed  was  increased. 

The  danger  arose  from  the  large  increment  in  the  aerodynamic 
forces  resulting  from  a  minor  pitch  angle  increment  at  high 
dynamic  pressures.  As  a  result,  the  model  tended  to  instantly 
fly  beyond  the  traveling  limits.  Without  an  active 
compensation  scheme  in  addition  to  the  manual  snubber  cables,  X 
the  consequence  could  be  catastrophic.  So  this  "kick"  was  not 

cs 

given  for  the  higher  airspeeds. 

For  comparison  purposes,  the  aeroelastic  behavior  of  the 
model  free  in  pitch  only  was  also  tested.  The  pitch  only 
configuration  was  set  up  by  blocking  the  model  plunging  freedom 
at  the  middle  of  the  test  section.  The  testing  variables  and 
procedure  were  the  same  as  those  in  the  testing  of  the  model  — 
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free  in  pitch  and  plunge.  Due  to  the  absence  of  the  plunge 
freedom,  the  risk  of  catastrophe  was  largely  reduced  in  this 
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CHAPTER  POUR 
RESULTS 


4.1  Model  Structural  and  Aerodynamic  Characteristics 


(a)  Basic  material  properties 


The  material  properties  of  Hercules  AS4/3501-6  prepreg  are 
summarized  in  Table  3.1.  The  longitudinal  flexural  modulus  E L 
was  obtained  from  the  results  of  the  vibration  test  of  the 
[02/90]s  wing.  The  longitudinal  flexural  modulus,  determined 
from  the  first  and  second  bending  frequencies  is  106  Gpa  which 
is  25%  lower  than  the  tensile  extensional  modulus  of  143  Gpa. 
This  stiffness  reduction  is  about  equal  to  that  of  ASl/3501-6 
reported  in  Ref.  40.  Since  no  testing  was  performed  to 
determine  the  actual  reduction  in  transverse  flexural  modulus, 
Et,  the  value  of  ASl/3501-6  prepreg  was  used.  This 
substitution  for  matrix  dominant  property  was  acceptable 
because  the  same  matrix  system  was  used  in  both  materials. 
Ref.  40  also  showed  that  the  shear  modulus  and  Poisson's  ratios 
are  almost  the  same  for  in-plane  and  out-of-plane  loadings. 
The  Hercules  AS4/3501-6  product  data  of  shear  modulus  and 
Poisson's  ratio  were  therefore  used  in  the  subsequent 
analysis.  For  comparison  purposes,  the  Hercules  AS4/3501-6 
product  data  of  in-plane  tensile  loading  are  also  summarized  in 
Table  3.1.  The  longitudinal  flexural  modulus  was  also 
determined  from  the  load-deflection  slope  of  the  cantilevered 
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[02/90]g  wing.  This  result  for  is  119  GPa, 


b)  Vibration  frequencies  and  modes 


Wing  alone  vibration  frequencies,  both  measured  and 
calculated,  are  given  in  Table  4.1.  The  frequency  variations 
with  respect  to  the  ply  orientation  agrees  with  the  previous 
investigation  in  Ref.  30.  For  unbalanced  laminated  wings,  the 
designations  of  "torsional  mode"  and  "bending  mode"  are 
somewhat  vague  because  of  stiffness  coupling. 

Figures  4.1  to  4.4  show  the  predicted  mode  shapes,  based 
on  the  seven-mode  Rayleigh-Ritz  analysis,  as  well  as  the 
observed  vibration  node  line  up  to  the  third  mode  for  all  wings 
tested.  For  the  (02/90]s  wing,  a  skewed  node  line  was  observed 
in  torsional  mode.  This  may  be  due  to  the  slight  stiffness 
coupling  because  of  the  misaligned  ply  orientation.  The  other 
possible  reason  for  the  skewed  node  line  is  due  to  the 
vibration  test  setup  for  a  30°  forward  swept  full-span  wing. 
Since  the  swept  wing  is  not  symmetric  with  respect  to  the 
excitation  point  (shaker  head),  the  inertia  force  contributions 
from  each  side  of  wing  can  affect  the  shaker  head  rigidity. 
Therefore,  the  observed  node  line  may  be  a  resultant  motion  of 
the  shaker  head  deformation  and  the  true  wing  torsional  mode. 
The  degree  of  abnormality,  however,  was  actually  amplified  by 
the  excessive  excitation  force  in  order  to  excite  the 
relatively  stiff  torsional  mode. 

To  assure  this  abnormality  remained  acceptable,  static 


Table  4.1  Vibration  Frequencies  of  Cantilever  Wings 


First 

First 

Second 

Bending 

Torsion 

Bending 

Wing 

calc . 

exp . 

calc . 

exp . 

calc . 

exp . 

(Hz) 

(Hz) 

(Hz) 

(HZ) 

(Hz) 

(Hz) 

[02/9°1s 

left 

right 

11.1 

11.4 

11.7 

36.2 

37.6 

37.1 

69 . 3 

71.8 

72.8 

tlV°U 

left 

right 

8.4 

9.0 

8.8 

46.7 

45.3 

44.9 

58.7 

64.3 

62.9 

1 302/0 ] g 

left 

right 

['152/01s 

5.6 

5.8 

6.0 

40.5 

37.4 

37.9 

60.6 

58.4 

57.7 

left 

right 


8.4  8.8 

8.9 
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46.5 
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62.6 
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Fig.  4.1  Cantilever  w  \  ng  Vibration  Modes,  [ 0 2/9 0 ) s  Wing 
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Fig.  4.2  Cantilever  wing  Vibration  Modes,  115-/0]  Wing 
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Fig.  4.4  Cantilever  wing  Vibration  Modes,  ( — 15-/0 )  wing 
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cantilever  test  was  performed  to  measure  the  bending  induced 
twisting  angle.  The  wing  was  bent  by  a  point  load  applied  at 
the  mid-chord  of  the  wing  tip  section.  The  wing  tip  deflection 
and  twisting  angle  were  then  mesaured  at  each  incremental 
loading  step.  Testing  results,  as  shown  in  Figure  4.5,  showed 
the  negligible  wing  tip  twisting  as  wing  bent.  This  test 
reassured  the  structural  integrity  of  the  ( 0^/90 ] s  wing. 

No  free-free  vibration  test  was  performed  on  the  aircraft 
model.  The  calculated  free-free  vibration  frequencies  and 
modes  are  tabulated  in  Table  4.2.  Comparison  between  the 
results  of  cantilever  and  free-free  wings  shows  the  slight 
increase  in  vibration  frequencies  for  the  free-free  results  as 
expected  due  to  the  inertia  coupling  of  rigid  body  modes  and 
wing  elastic  modes.  All  the  vibration  analyses  include  the 
virtual  mass  effect  so  that  the  results  are  comparable  with 
actual  vibration  tests. 
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TIP  TWISTING  (deg) 


TIP  DEFLECTION  (xlO“2  m) 


Fig.  4.5  Cantilever  Load  Deflection  Curve,  {02/90)s  Wing 
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Table  4.2  Calculated  Fundamental  Frequency  and 
Vibration  Mode  of  Free-Free  Wings 


[02/90]s 

[152/0]s 

[302/0]s 

[ -152/0 ] 

eq. 

11.3  Hz. 

8.6  Hz. 

6.0  Hz. 

8.6  Hz 

*1 

-.04 

-.04 

-.04 

-.04 

<*2 

.13 

.13 

.12 

.13 

^3 

1.0 

1.0 

1.0 

1.0 

<*4 

.00 

-.62 

-.44 

.62 

.00 

-.01 

-.01 

-.01 

^6 

.00 

-.05 

-.04 

.05 

^7 

-.02 

-.04 

-.04 

-.04 

q^:  plunge  mode 
q2 :  pitch  mode 
q^:  first  bending  mode 
q4:  first  torsional  mode 
q^:  second  bending  mode 
qg :  second  torsional  mode 
q^:  chordwise  bending  mode 


rig.  4.6  Lift  Coefficient  for  Aircraft  Model 
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lift  coefficient  plot,  the  canard  contributions  to  the  overall 
lift  are  not  really  noticeable.  The  measurements  of 
aerodynamic  moment  coefficients,  however,  experienced  greater 
variations  in  the  linear  portion  and  stall  regions.  These 
variations  in  the  linear  range  were  possibly  caused  by  the 
relatively  higher  pitch  friction  incurred  in  the  measuring 
procedure,  wing  warping,  and  gage  wire  attachments  leading 
outside  the  fuselage.  The  extreme  sensitivity  to  the  wing 
warping  of  aerodynamic  moment  is  also  evident  in  the  results  of 
varied  canard  setting.  As  shown  in  Figure  4.7,  the  aerodynamic 
moment  coefficients  are  significantly  changed  as  canard  setting 
was  moderately  varied.  The  variation  in  the  stall  region  was 
largely  due  to  the  nonlinearity  of  aerodynamic  interference 
effect . 

Based  on  Figures  4.6  and  4.7,  the  results  show  values  of 

CL  -  5.33  and  CM  •  -0.97  which  give  the  present  aircraft  model 
a  a 

an  18%  static  margin,  calculated  from  -CM  /CL  .  In  addition, 

a  a 

the  model  pitch  decay  rates  in  still  air  and  at  5  m/sec  were 

calculated  to  give  the  CM>.  The  result  shows  value  of  CM." 

a  a 

-4.95. 


4.2  Cantilever  Wing  Divergence  and  Flutter 


For  comparison  purposes,  the  flutter  and  divergence 
boundaries  of  the  cantilever  wing  were  tested  and  are  presented 
in  Figure  4.8.  Because  of  the  wing  flexibility,  the 
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Fig.  4.8  Cantilever  Wing  Flutter  and  Divergence  Boundaries 
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aerodynamic  nonlinearity,  i.e.,  wing  stall,  would  occur  before 
the  wing  model  was  overstressed  or  failed.  Aeroelastic 
instability  was  thus  encountered  and  observed  in  the  form  of 
limit  cycle  oscillation.  Linear  aeroelastic  phenomena,  i.e., 
wing  divergence  and  bending-torsion  wing  flutter,  were 
identified  at  low  fuselage  (wing  root)  angles  of  attack. 

The  cantilever  wing  results  shown  in  Figure  4.8  are 
similar  to  those  given  in  Ref.  10.  For  the  (152/0]s  and 
(302/0]s  wings,  bending-torsion  flutter  occured  at  low  fuselage 
angles  of  attack.  At  higher  angles  of  attack,  the 
bending-torsion  wing  flutter  gradually  changed  to 
torsional-stall  flutter  (TSF)  with  an  accompanying  increase  in 
the  flutter  frequency  toward  the  still  air  torsional  mode 
frequency,  and  a  decrease  in  the  flutter  speed.  Calculations 
of  the  linear  bending-torsion  flutter  speed  and  frequency 
showed  good  agreement  with  experiment  at  low  angles  of  attack. 
In  the  case  of  the  [ 1 52/O 1 s  wing,  the  bending-tor sion  flutter 
speed  at  zero  angle  of  attack  was  30  m/sec  compared  with  the  28 
m/sec  calculated  flutter  speed.  The  measured  and  calculated 
flutter  frequencies  were  31  and  30  Hz  respectively.  As  the 
fuselage  angle  of  attack  increased  from  08  to  16°,  a  46%  drop 
in  flutter  speed  was  observed.  In  the  case  of  the  the  [ 30^/0  ]  s 
wing,  at  tunnel  speed  limit  of  30  m/sec,  flutter  occured  only 
when  fuselage  angle  of  attack  was  beyond  7.5°.  No  flutter  was 
observed  with  a  smaller  angle  of  attack  at  this  speed  limit. 
The  predicted  bending-torsion  flutter  speed  is  38  m/sec  which 
is  beyond  the  tunnel  speed  limit.  The  measured  torsional-stall 
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flutter  speed  dropped  40%  as  the  wing  root  angle  of  attack 
increased  from  6°  to  15®. 

For  the  (02/901s  and  [-152/0Js  wings,  divergence  occured 
at  low  fuselage  angles  of  attack.  This  divergence  condition 
was  noted  when  the  fuselage  angle  of  attack  could  not  be  set 
small  enough  to  keep  the  wing  from  "flipping  over"  to  either 
side.  As  the  fuselage  angle  of  attack  was  increased,  the 

static  divergence  quickly  changed  to  either  a  torsional-stall 
flutter  ([02/90Js  wing)  or  a  bending-stall  flutter  ( BSF ) 

([-152/0]  wing),  with  an  accompanying  decrease  in  flutter 
speed.  In  the  case  of  the  [ 0 2/9 0 ] s  wing,  the  calculated 
divergence  speed  is  21  m/sec  compared  with  the  21  m/sec 
measured  wing  divergence  speed.  When  the  divergence  point  was 
encountered,  however,  the  [02/901s  win<3  still  exhibited  a  45  H2 
torsional  oscillation  in  addition  to  the  phenomenon  of 
"flipping  over".  This  concurrence  is  due  to  the  fact  that 
divergence  speed  is  only  slightly  lower  than  the  flutter  speed 
predicted  at  23  m/sec.  Like  the  previously  mentioned  flutter 
prone  wings,  the  torsional-stall  flutter  existed  at  higher 
angles  of  attack  and  continued  its  presence  as  the  "static” 
divergence  was  encountered  at  low  angle  of  attack.  The 

observed  flutter  speed  dropped  52%  as  the  fuselage  angle  of 
attack  was  changed  from  0°  to  15°.  In  the  case  of  the 

[  —  1 5 2 0  )  s  wing,  the  divergence  speed  was  noted  at  14  m/s ec 
compared  with  the  13  m/sec  calculated  result.  It  should  be 
noted  that  the  so-called  static  divergence  actually  appeared  as 
a  dynamic  bending-stall  flutter.  As  fuselage  angle  of  attack 


was  increased,  the  bending-stall  flutter  speed  slowly  decreased 
while  its  flutter  frequency  slowly  increased  toward  the  still 
air  wing  bending  frequency.  Since  the  phenomenon  of 
bending-stall  flutter  persisted  from  low  to  high  angle  of 
attack,  the  (-152/0)s  wing  exhibits  the  most  smooth  transition 
with  angle  of  attack  change.  In  comparison,  the  flutter  prone 
wings,  ( 1 52/O ] s  and  { 3 0 0 ) s  showed  a  more  dramatic  speed  drop 
due  to  the  transition  from  one  phenomenon,  bending-torsion 
flutter,  to  the  other,  torsional-stall  flutter. 


4 . 3  Body  Freedom  Flutter  ( BFF )  and  Support  Instability 

Body  freedom  flutter  was  predicted  by  the  linear  analysis 
and  was  observed  in  the  wind  tunnel  tests.  There  also  existed, 
however,  an  additional  wind  tunnel  support  related  dynamic 
instability  involving  the  model  plunge  and  pitch  motions. 
Since  the  support  instability  occurred  before  the  body  freedom 
flutter,  it  was  of  concern  to  the  integrity  of  the  aeroelastic 
behavior  and  the  overall  testing  objective.  Therefore,  before 
presenting  the  body  freedom  flutter  results,  this  support 
instability  will  first  be  described. 

( a )  Support  instability 


The  aircraft  model  on  the  original  support  system  had  a 
0.63  Hz  plunging  frequency  which  resulted  from  the  suspension 
spring,  but  had  no  mechanical  pitching  stiffness  except  for  a 
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small  amount  of  pendulum  stiffness  which  gave  a  0.20  Hz 
pitching  frequency  in  still  air.  Measurements  from  aircraft 
model  step  responses  gave  values  of  support  plunging  and 
pitching  equivalent  viscous  damping  ratios  of  approximately 
0.06  for  plunge  motion  and  0.26  for  pitch  motion.  With  the 
averaged  damping  ratios,  a  simplified,  two  degrees  of  freedom 
longitudinal  stability  analysis,  discussed  in  Sec.  2.3,  showed 
that  a  support  instability  would  occur  at  8  m/sec  and  end  at  31 
m/sec.  A  useful  tool  showing  the  interaction  of  plunging  and 
pitching  motions  is  to  plot  the  frequency  components  against 
the  airspeed,  namely  an  oo-V  plot.  Figure  2.7  shows  the 
predicted  interaction  for  a  "rigid"  wing  aircraft  model. 

The  complete  7-mode  aeroelastic  analyses,  as  discussed  in 
Sec.  2.2,  showed  that  wing  flexibility  could  significantly 
modify  the  termination  speed  of  the  support  instability,  but 
not  the  onset  speed.  Table  4.3  shows  the  predicted  and 
observed  support  instability  boundaries  and  their  good 
agreement  for  all  four  wings.  The  frequencies  of  the  support 
instability  were  in  the  range  of  0.5  to  0.6  Hz  which  is  near 
the  vibration  frequency  of  the  support  plunge  motion  in  still 
air.  Figure  4.9  shows  an  example  of  support  instability  -'te:* 
the  model  was  seen  to  plunge  unboundedly  until  it  mt  *“e  * 
and  bottom  plunge  stops,  and  the  wing  was  see"  *  •  '■ » 

"rigid"  until  high  model  pitch  angle  was 
records  of  support  instability  versus  a  .  .  ■  > 

also  shown  in  Figure  4.10. 

When  the  support  i  r.  s  t  a  r  i  1  .  '  .  -  a  -  • 


Table  4.3  Support  Instability  Boundaries 
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(302/0]s 


calc .  exp . 


m/sec 

(Hz) 

m/sec 
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two  factors  vital  to  the  survivability  of  the  ongoing  testing. 
One  was  the  time  constant  (time  to  double  the  amplitude)  of  the 
unstable  mode,  the  other  was  the  model  robustness.  Without  an 
active  compensation  scheme,  an  instability  with  a  short  time 
constant  could  not  be  limited  by  manually  operated  snubber 
cables.  The  unstable  mode  exhibited  in  the  original  support 
system,  fortunately,  appeared  to  be  mild  and  with  a  long  time 
constant,  and  therefore  could  be  guarded  and  flown  through  if 
the  model  was  not  given  much  of  an  initial  disturbance.  This 
additional  nonlinear  phenomenon  of  disturbance  magnitude 
dependency,  which  arose  from  friction  in  the  support  system, 
was  often  observed  in  the  higher  airspeed  range  and 
occasionally  in  the  case  of  using  a  linear  bearing  which  was 
contaminated  by  particles  amid  wind  tunnel  airflow  after  long 
hours  of  running. 

To  suppress  the  wind  tunnel  support  related  dynamic 
instability,  the  original  pitching  support  was  modified  by 
adding  mechanical  pitching  springs,  so  that  the  pitching 
frequency  was  increased  from  0.20  Hz  to  0.85  Hz  as  compared  to 
the  plunging  frequency  of  0.63  Hz.  The  modified  support 
pitching  damping  ratio  was  now  reduced  to  0.03  due  to  the 
increased  pitching  frequency.  This  modified  support  system  was 
tested  and  showed  no  signs  of  the  support  instability 
previously  encountered. 


The  primary  interest  in  flutter  research  usually  centers 
on  the  flutter  boundary  and  frequency  of  oscillation,  while  the 
flutter  mode  and  subcritical  behavior  are  often  seen  as  of 
secondary  interest.  In  the  present  study,  the  flutter  mode  and 
transient  behavior  were  both  monitored  and  played  important 
roles  in  verifying  the  instability  mechanism  and  assessing  the 
preliminary  aircraft  flying  quality.  In  this  section,  the  body 
freedom  flutter  testing  results  will  be  presented  in  two 
stages.  First,  the  general  behavior  will  be  described.  Then, 
the  results  of  each  model  will  be  detailed  in  following. 

Body  freedom  flutter  was  observed  during  the  tests 
wherever  predicted  analytically.  The  identification  of  the 
body  freedom  flutter  boundary  depended  on  the  observation  of 
model  transient  behavior.  First,  visual  observation  was  made 
during  the  test.  The  oscillograph  records  of  rigid  body  and 
wing  responses  were  then  examined  in  terms  of  frequency  and 
phase  relations.  Lastly,  video  movies  of  the  chronological 
events  were  replayed  for  flutter  confirmation.  Figure  4.11 
shows  a  typical  flutter  onset  of  :he  model  with  the  l 0 2^^ ®  1  s 
wing.  In  this  figure,  the  aeroelastic  responses,  prior  to  and 
on  the  flutter  onset,  were  seen  to  be  in  sharp  contrast  for 
this  particular  wing.  The  strength  of  instability,  however, 
can  be  better  understood  via  the  results  of  root  locus  analysis 
which  will  be  discussed  later.  Even  though  there  were  pitch 
decay  rates  estimated  from  subcritical  measurements,  the 
results  showed  too  much  scattering  and  insufficient  convergence 
as  shown  in  Figure  4.12,  and  could  not  be  used  to  correlate  the 
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instability  onset.  It  is  interesting  to  note  that  the  results 
of  model  free  in  pitch  only  actually  are  fairly  good  and  can  be 
used  to  correlate  the  flutter  results. 

During  the  test,  the  body  freedom  flutter  was  observed  to 
be  mild  and  in  the  form  of  a  limit  cycle  oscillation.  The  mild 
instability  onset  was  partially  a  reflection  of  the  low 
frequency  (1.0  Hz  to  2.8  Hz)  of  the  instability.  The  limiting 
amplitude  of  the  limit  cycle  oscillation  was  undoubtedly  due  to 
the  aerodynamic  nonlinear  stalling  phenomenon,  rather  than  the 
structural  nonlinearity.  These  features  proved  to  be  useful 
since  it  made  close  examination  of  the  flutter  modes  and 
nonlinear  behavior  possible. 

All  the  experimental  flutter  modes  exhibited  a  nodal  point 
near  the  model's  nose  position,  while  wing  bending  was 
approximatly  in-phase  with  the  model  pitch  motion,  and  the 
plunge  motion  was  about  180°  out-of-phase.  Wing  torsional 
response  was  dictated  by  the  inherent  elastic  coupling  of  the 
composite  wings.  The  1 0 2/90 ] s  wing  exhibited  the  typical 
aeroelastic  effect  of  a  forward  swept  wing  with  the  wing's 
upward  bending  accompanied  by  leading  edge  up  twisting.  In  the 
case  of  wing,  the  leading  edge  up  twisting  was 
amplified  by  the  adverse  elastic  coupling.  For  the  [ 152/0 ] s 
and  I  3 0 2/O ] s  wings,  the  elastic  coupling  worked  in  an  opposite 
way  such  that  the  torsional  gages  showed  a  leading  edge  down 
twisting  as  the  wing  bent  upward.  Table  4.4  summarizes  the 
corresponding  body  freedom  flutter  mode  predictions  in  terms  of 
the  contributions  of  the  uncoupled  assumed  modes. 


-104- 


Table  4.4  Calculated  Body  Freedom  Flutter  Modes 

(with  modified  support  system) 


Wing 

(o2/9°]s 

(302/0]s 

( -152/0 ] 

*1 

.  45/-17  5 0  * 

. 68/-171 0 

. 42/-174 

^2 

1.0/  0° 

1.0/  0° 

1.0/  0“ 

<*3 

. 27/-13  8 

. 87/-2  4  ° 

. 22/-1 4  0 

<*4 

. 12/-12  0 

.34/154° 

. 16/-1 4  ° 

^5 

0 

.01/21° 

0 

^6 

0 

.05/159° 

.01/-110 

^7 

0 

.04/154“ 

.01/166“ 

q.^ :  plunge  mode 
q2 :  pitch  mode 
q^:  first  bending  mode 
q^:  first  torsional  mode 
q^:  second  bending  mode 
q^ :  second  torsional  mode 
q-y :  chordwise  bending  mode 

* 


amplitude/phase  angle 


I 


In  the  case  of  the  model  free  in  pitch  only,  the  flutter 
modes  were  similar  to  those  previously  mentioned  except  that 
the  model  was  forced  to  pitch  about  a  fixed  pitching  axle. 
Table  4.5  summarizes  the  calculated  body  freedom  flutter 
modes . 

Another  concern  in  observing  the  flutter  boundary  was  the 
resemblance  of  the  flutter  mode  and  gust  response.  In  the  case 
of  gust  response,  the  wing  will  bend  up  elastically  for  a  given 
nose-up  pitch  angle.  Therefore,  the  behavior  recorded  at  the 
body  freedom  flutter  onset  would  be  similar  to  the  turbulent 
gust  response  or  the  intentional  "kick"  at  subcritical  speed, 
when  the  body  freedom  flutter  onset  was  approached,  the  gust 
response  became  so  lightly  damped  that  the  determination  of 
body  freedom  flutter  must  be  repeatedly  examined.  As  a  result, 
the  flutter  point  was  formally  determined  only  after  a  harmonic 
oscillation  was  fully  developed. 

In  the  following,  flutter  behavior  of  all  models  tested 
will  be  described  in  detail.  For  the  model  with  the  [ 0 2/^ 0 ] s 
wing,  the  experimental  body  freedom  flutter  mode  was  usually 
accompanied  by  a  small  high  frequency  torsional  oscillation  as 
shown  in  Figure  4.13.  This  torsional  oscillation  was  a 
localized  phenomenon  and  generally  occurred  at  the  peaks  of 
body  freedom  flutter  oscillation.  This  torsional  frequency  was 
close  to  that  of  the  cantilever  wing  torsional-stall  flutter 
shown  before  and  occured  as  the  model  went  into  the  higher 
pitch  angles  of  the  body  freedom  flutter  oscillation.  Since 


the  period  of  the  body  freedom  flutter  was  much  longer  than 


Table  4.5  Calculated  Body  Preedoa  Flutter  Modes 

(model  free  in  pitch  only) 


Wing 

[02/9°)s 

(302/°)s 

( -152/0 ] 

<*2 

1.0/  0°* 

1.0/  0° 

1.0/  0° 

<*3 

. 18/-12  0 

. 52/-16  0 

. 15/-110 

** 
t T 

. 08/-1 3  0 

.20/163° 

.11/-110 

<*5 

0 

0 

0 

0 

.03/165“ 

.  01/-9  ° 

<*7 

0 

.02/162° 

0 

q2 :  pitch  mode 
q^:  first  bending  mode 
q4:  first  torsional  mode 
q^:  second  bending  mode 
qg :  second  torsional  mode 
q^ :  chordwise  bending  mode 

* 


amplitude/phase  angle 


bending 

strain 

(xlO-3) 


+  2 


torsion 

strain 


that  of  the  torsional-stall  flutter,  the  latter  could  build  up 
several  cycles  before  the  model  pitch  angle  became  small 
again.  The  body  freedom  flutter  speed  and  flutter  frequency 
are  summarized  in  Table  4.6  for  all  configurations  tested, 
i.e.,  model  in  original  and  modified  support  systems,  and  model 
free  in  pitch  only. 

For  the  model  with  the  wing,  no  body  freedom 
flutter  was  encountered  within  the  30  ra/sec  tunnel  speed  limit, 
except  for  the  model  free  in  pitch  only  which  exhibited  a  1.4 
Hz  body  freedom  flutter  at  22  m/sec.  In  steadily  trimmed 
flight,  without  the  presence  of  body  freedom  flutter,  the  wing 
alone  exhibited  a  bending-torsion  flutter  at  29  m/sec,  as  shown 
in  Figure  4.14,  which  became  the  critical  unstable  mode. 
Comparing  with  the  cantilever  wing  results  shows  that  the 
effect  of  rigid  body  freedoms  on  bending-torsion  wing  flutter 
are  not  noticeable.  Table  4.7  summarizes  the  body  freedom 
flutter  testing  results. 

For  the  model  with  the  l^O^/OJg  wing,  the  observed  body 
freedom  flutter  mode,  as  shown  in  Figure  4.15,  was  similar  to 
that  of  the  ^2/90)s  wing  except  that  no  torsional-stall 
flutter  appeared  in  the  body  freedom  flutter  oscillation.  This 
result  was  consistent  with  the  cantilever  wing  testing  in  which 
no  torsional-stall  flutter  occured  within  the  tunnel  speed 
limit,  unless  a  high  pitch  angle  was  present.  Table  4.8 
summarizes  the  body  freedom  flutter  speed  and  flutter  frequency 
for  all  the  configurations  tested. 

For  the  model  with  the  (-15, /0]  wing,  a  clean 
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Table  4.6  Body  Freedoa  Flutter  Boundary,  (02/90]s  Wing 


Model  on  Original  Support 


calc . 

exp. 

*c 

run# 

m/sec 

(Hz) 

m/sec 

(Hz  ) 

0.0° 

185 

19. 

(2.7) 

19. 

(2.6) 

2.5° 

180 

19. 

(2.7) 

5.0° 

181 

19. 

(2.6) 

Model  on  Modified  Support 

calc . 

exp . 

&c 

run# 

m/sec 

(Hz) 

m/sec 

(Hz) 

2.5° 

127 

19. 

(2.7) 

20. 

(2.8) 

Model  Free  in 

Pitch  Only 

calc. 

exp. 

*c 

run# 

m/sec 

(Hz) 

m/sec 

(Hz) 

0.0° 

168 

17. 

(2.0) 

16. 

(1.7) 

2.5° 

169 

17. 

(1.9) 

5c  : 

canard 

setting 

A 
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Table  4.7  Body  Freedom  Flutter  Boundary,  (152/018  wing 


Model  on  Original  Support 

calc . 

exp. 

*c 

run# 

m/sec  (Hz) 

m/sec  (Hz) 

2.5° 

198 

* 

— 

5.0° 

199 

* 

Model  on  Modified  Support 

calc . 

exp. 

*c 

run# 

m/sec  (Hz) 

m/sec  (Hz) 

2.5° 

133 

* 

- 

Model  Free  in 

Pitch  Only 

calc . 

exp. 

Sc 

run# 

m/sec  (Hz) 

m/sec  (Hz) 

0.0° 

110 

* 

22.  (1.4) 

2.5° 

112 

★ 

22.  (1.4) 

2.5° 

113 

* 

20.  (1.5) 

Sc  s 

canard 

setting 

*  :  No  BFF  predicted  analytically 

:  Wing  flutter  was  encountered 
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JC. 
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.V 
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Table  4.8  Body  Freedoa  Flutter  Boundary,  HOj/Olg  Wing 

Model  on  Original  Support 


calc . 

exp . 

Sc 

run# 

m/sec  ( Hz ) 

m/sec 

(Hz) 

0.0° 

207 

28.  (2.6) 

n .  a . 

2.5° 

206 

26. 

(1.7) 

5.0° 

211 

n .  a . 

Model  on  Modified  Support 

calc . 

exp . 

*c 

run# 

m/sec  (Hz) 

m/sec 

(Hz) 

2.5° 

137 

28.  (2.6) 

27. 

(1.8) 

Model  Free  in 

Pitch  Only 

calc . 

exp . 

Sc 

run# 

m/sec  (Hz) 

m/sec 

(Hz) 

0.0° 

129 

21.  (1.8) 

19. 

(1.2) 

2.5° 

130 

19. 

(1.2) 

2.5° 

133 

19. 

(1.3) 

6c  5 

canard 

setting 

n .  a . 

:  data 

not  available 

a 


«•, 

j*. 


i 

V 


.v 


3 

.v 


pitch/bending  oscillation  of  body  freedom  flutter  was  observed 
at  the  flutter  onset  as  shown  in  Figure  4.16.  The  body  freedom 
flutter  speed  and  frequency  are  tabulated  in  Table  4.9. 

As  previously  mentioned  in  cantilever  wing  testing,  the 
static  divergence  of  the  (— 1S2/01S  wing  was  actually  observed 
as  a  dynamic  bending-stall  flutter  with  a  somewhat  lowered 
bending  frequency.  When  the  wing  root  boundary  condition  was 
relaxed,  it  was  important  to  be  able  to  distinguish  the  body 
freedom  flutter  from  the  bending-stall  flutter.  From  the 
observations  of  the  model's  behavior,  there  were  three  features 
which  existed  for  distinction  between  the  two  instabilities. 
First,  the  observed  body  freedom  flutter  mode  was  similar  to 
that  of  models  with  the  [302/0]s  and  the  [02/90)s  wings  which 
did  not  experience  any  bending-stall  flutter  in  the  cantilever 
wing  testing. 

Secondly,  the  body  freedom  flutter  oscillated  in  a 
relatively  low  frequency,  2.3  Hz  for  the  model  with  the 
[ — 1 5 2/0 ) s  wing,  compared  to  the  5.0  Hz  bending-stall  flutter 
observed  in  the  cantilever  wing  testing.  Lastly,  as  the  tunnel 
airspeed  was  increased  beyond  the  body  freedom  flutter  onset, 
the  wing  bending  response  was  increased  and  accompanied  by 
trace  of  bending-stall  flutter  at  each  peaks  of  the  well 
developed  body  freedom  flutter  oscillation.  The  frequency  of 
this  localized  bending-stall  flutter  was  found  to  be  close  to 
the  frequency  of  the  cantilever  wing  bending-stall  flutter. 
Therefore,  the  bending-stall  flutter  was  concluded  not  to 
appear  on  a  "free  flight"  testing  except  for  the  localized 
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Table  4.9  Body  Freedoa  Flutter  Boundary,  ( -152/0  J  g  Wing 

Model  on  Original  SuDDort 


calc 

exp. 

m/sec 

(Hz) 

m/sec 

( 

12. 

(1.8) 

13. 

(2 

13. 

(2 

Model 

on  Modified  Support 

calc 

exp . 

m/sec 

( Hz ) 

m/sec 

( 

12. 

(1.9) 

13. 

( 

Model 

Free  in 

Pitch  Onl 

canard  setting 


appearance  in  the  body  freedom  flutter  oscillation. 

Unlike  the  torsional-stall  flutter,  the  bending-stall 
flutter  intensified  quickly  with  a  small  increase  in  airspeed. 
The  severe  wing  flapping  motion  in  bending-stall  flutter  could 
induce  significant  inertia  force  in  plunge  motion  to  detach  the 
friction-fitted  linear  bearing  from  the  model,  and  interrupt 
the  ongoing  testing.  In  some  of  these  cases,  when  model  roll 
freedom  was  relaxed,  due  to  loss  of  the  bearing,  and  when  the 
pitch  angle  was  held  by  snubber  cables,  the  [—15 2/0 ] s  wing 
could  exhibit  either  a  6  Hz  symmetric  bending-stall  flutter  or 
an  11  Hz  anti-symmetric  bending-stall  flutter. 


CHAPTER  FIVE 


DISCUSSION 


In  this  chapter,  the  significant  trends  of  the  previous 
results  will  be  summarized  first.  Then  a  comparison  will  be 
made  with  theoretical  analysis  for  the  flutter  boundary  and 
transient  aeroelastic  responses.  With  this  well-verified 
analytical  model,  some  parametric  effects  on  body  freedom 
flutter  will  be  explored.  Lastly,  the  aeroelastic  effect  on 
support  instability  will  be  further  examined. 


5 . 1  Summary  of  Results 

Typical  results  of  the  various  wings  tested  are  summarized 
in  Table  5.1.  Figure  5.1  also  gives  an  illustrative  summary  of 
the  flutter  boundary  in  terms  of  aeroelastic  tailoring.  From 
these  results,  four  major  trends  are  noted  in  following. 

( a )  Aeroelastic  tailoring  effect 

The  aircraft  model  with  the  [15 2/® 1 s  w*-n9  demonstrated  the 
best  tailored  aeroelastic  behavior  for  all  wings  tested.  From 
Figure  5.1,  it  can  be  seen  that  the  "optimum"  tailored  wing,  in 
terms  of  body  freedom  flutter,  has  a  ply  orientation  between 
15°  and  20°.  Due  to  the  presence  of  the  other  type  of 
instability,  namely,  wing  bending-torsion  flutter,  the  actual 


Table  5.1  Suaaary  of  Divergence  and  Flutter  Boundaries 


i 

■ 

» 

b 

i. 


Wing  (  02/90Js  ( 1 5  2/0  J  s 

calc.  exp.  calc.  exp. 

m/sec  (Hz)  m/sec  (Hz)  m/sec  (Hz)  m/sec  (Hz) 

Cantilever  j 

Wing  21  21  ***  - 

Divergence 

Cantilever 

Wing  23  (23.)  ---  28  (30.)  30  (31.) 

Flutter 


( original 
support ) 

19 

(2.7) 

19 

(2.7) 

★  ★  * 

>30 

BFF 

( modified 
support ) 

19 

(2.7) 

20 

(2.8) 

★  *  * 

>30 

BFF 

(free  in 
pitch  only) 

BFF 

17 

(2.0) 

17 

( 1 .  V  ) 

} 

★  *  * 

22 

( completely 
unrestrained ) 

19 

(2.6) 

n .  a 

★  *  * 

n . 

1  No  BFF  predicted  analytically 

2  Wind  tunnel  speed  limit  -  30  m/sec 

3  Not  applicable 


Table  5.1  Suoaary  of  Divergence  and  Flutter  Boundaries 

( cont ' d ) 


Wing  (  302/0  Js  [—15  2/° 1 s 

calc.  exp.  calc.  exp. 

m/sec  (Hz)  m/sec  (Hz)  m/sec  (Hz)  m/sec  (Hz) 

Canti leve  r 

Wing  45  -  13  14 

Divergence 

Cantilever  , 

Wing  38  (  3  3  .  )  >30  ***  - 

Flutter 


BFF 

(original  28  (2.6)  26  (1.7)  12  (1.8)  13  (2.3) 

support ) 

BFF 

(modified  28  (2.6)  27  (1.8)  12  (1.9)  13  (2.3) 

support ) 

BFF 

(free  in  21  (1.8)  19  (1.2)  11  (1.4)  10  (1.1) 

pitch  only) 

BFF  3 

(completely  27  (2.5)  n.a.  12  (1.7)  n.a. 

unrestrained ) 


1  No  BFF  predicted  analytically 

2  Wind  tunnel  speed  limit  -  30  m/sec 

3  Not  applicable 


VELOCITY  (m/sec) 


calc . 


Cantilever  Wing  Flutter 
Cantilever  Wing  Divergence 
Body  Freedom  Flutter 


unstable 


/V 

unstable  y  y 


/ 


stable 


^unstable 

\ 


PLY  ANGLE  (deg) 


'02/0' 
»  « 


Body  Freedom  Flutter  Boundary, 
Model  on  Modified  Support  System 
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"optimum"  ply  orientation  may  be  altered.  Examination  of 
Figure  5.1  shows  that  the  wing  flutter  speed  can  also  be 
improved  by  rotating  the  ply  orientation.  However,  with  the 
same  ply  rotation  angle,  the  improvement  of  wing  flutter  seems 
relatively  mild  compared  to  the  body  freedom  flutter.  In  the 
present  study,  the  overall  "optimum"  ply  orientation  is  only 
slightly  changed,  and  is  in  the  vicinity  of  25®.  As  expected, 
the  model  with  the  [— 152/01s  wing  has  the  most  deteriorated 
performance  because  of  the  adverse  elastic  coupling  as 
discussed  in  Chapter  2.  Even  though  the  [  3 0 2/0  )  s  wing  has  the 
highest  coupling  stiffness  ratio,  D^g//  bgg  of  0.80,  the 
body  freedom  flutter  speed  actually  drops  instead  of 
increases.  This  is  due  to  the  fact  that  the  wing  bending 
stiffness  ratio,  Dn/D66'  simutaneously  decreased  as  the 
coupling  stiffness  ratio  increases  for  this  particular 
laminate.  The  adverse  effect  of  reduced  bending  stiffness 
overwhelms  the  improving  effect  of  stiffness  coupling. 


( b )  Body  freedom  flutter  vs.  wing  divergence 


Comparing  with  the  cantilever  wing  divergence  boundary, 
body  freedom  flutter  is  shown  to  be  more  critical  than  wing 
divergence  in  Figure  5.1.  The  degree  of  criticality,  however, 
is  not  uniform  for  all  tailored  wings.  In  other  words,  the 
aeroelastic  tailoring  effects  on  wing  divergence  and  body 
freedom  flutter  are  somewhat  different.  For  example,  the 
calculated  wing  divergence  speed  increased  by  114%  as  ply 
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orientation  rotates  from  0°  to  30°,  while  the  body  freedom 
flutter  speed  improved  by  only  421.  In  general,  the  "optimum" 
tailored  wing  seems  to  be  more  sensitive  to  the  presence  of 
rigid  body  freedoms. 

Examination  of  the  body  freedom  flutter  boundary  in 
Figure  5.1  shows  that  the  effect  of  aeroelastic  tailoring  on 
body  freedom  flutter  follows  the  trend  of  cantilever  wing 
divergence  tailoring.  This  suggests  that  wing  destiffening  is 
still  the  primary  triggering  mechanism  of  the  body  freedom 
flutter  as  of  the  cantilever  wing  divergence.  The  presence  of 
rigid  body  freedoms  modify  the  boundary  conditions  of  the 
aeroelastically  destiffened  wing  mode,  and  ultimately  changes 
the  instability  nature  from  the  static  divergence  to  dynamic 
flutter.  The  recognition  of  the  "destiffened"  wing  mode 
(analogy  to  wing  divergence  mode)  in  body  freedom  flutter  is 
not  only  important  in  understanding  the  underlying  instability 
mechanism  but  also  vital  to  the  possible  instability 
preventation . 

The  significance  of  the  above  observation  is  that  the 
prevention  of  body  freedom  flutter  will  largely  rely  on  the 
control  of  the  "destiffened"  wing  mode  instead  of  the  flight 
dynamic  modes,  even  though  these  involve  significant  rigid  body 
motions  in  body  freedom  flutter.  The  aeroelastic  tailoring  in 
the  present  study  serves  exactly  this  purpose  through  the  use 
of  composite  materials.  It  is  interesting  to  note  that  Ref.  41 
showed  weak  performance  of  the  flight  dynamics  control  by  using 
a  canard  based  stability  augmentation  system.  Refs.  42  and  43 
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also  demonstrated  the  possible  application  of  wing  based  active 
control  of  the  "desti f fened"  wing  mode.  Since  the  aeroelastic 
destiffening  is  a  strong  instability  tendency  while  the  short 
period  mode  is  originally  stable,  the  dominance  of  the 
"desti f fened"  wing  mode  in  body  freedom  flutter  would  be 
expected . 

Qualitatively  speaking,  the  aeroelastic  tailoring  trends 
of  the  cantilever  wing  divergence  can  be  used  as  a  preliminary 
design  for  preventing  the  body  freedom  flutter.  However,  the 
tailoring  effectiveness  is  not  quite  the  same  for  these  two 
types  of  instability.  The  actual  body  freedom  flutter 
boundary,  which  is  the  most  critical  now,  can  only  be  obtained 
by  the  complete  aeroelastic  analysis  including  the  rigid  body 
freedoms . 


(c)  Pitch-plunge  testing  vs.  pitch  only  testinc 


The  absence  of  the  plunging  freedom  resulted  in  a  somewhat 
lowered  body  freedom  flutter  speed  in  the  present  study, 
thereby  giving  a  conservative  flutter  boundary  as  shown  in 
Figure  5.2.  The  stabilizing  effect  of  plunging  freedom  is 
similar  to  the  findings  in  Ref.  20.  A  possible  explanation  is 
that  the  model  is  stabilized  by  the  favorable  phase  response  of 
the  pitch  and  plunge  freedoms.  When  the  aircraft  was  free  to 
pitch  and  plunge,  the  flutter  mode  always  exhibits  a  nodal 
point  near  the  model  nose  position  which  is  far  ahead  of 


pitching  axle.  The  longer  distance  between  the  nodal  point  and 
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airccaft  aerodynamic  center,  slightly  aft  of  the  pitching  axle, 
may  be  the  major  reason  for  delaying  the  instability.  The 
resulting  flutter  frequencies,  ranging  from  1.4  to  2.0  Hz,  were 
also  found  to  be  significantly  lower  than  those  of  the  models 
free  in  pitch  and  plunge  freedoms. 

Wind  tunnel  testing  of  the  model  free  in  pitch  only  is  a 
relatively  easy  task  and  an  attractive  alternative  for  body 
freedom  flutter  testing.  The  aforementioned  comparisons  of 
flutter  boundary  and  flutter  behavior,  however,  show  that  the 
aircraft  plunge  freedom  plays  an  important  role  in  body  freedom 
flutter  instability.  To  obtain  meaningful  results  and  true 
insight  of  body  freedom  flutter  of  a  forward  swept  wing 
aircraft,  a  more  complicated  testing,  including  both  rigid  body 
pitch  and  plunge  freedoms,  must  be  performed  such  as  in  the 
present  study. 

( d )  Analytical  Correlation 

The  aeroelastic  analysis  of  the  wind  tunnel  model  was 
based  on  the  analytical  model  formulated  in  Chapter  2.  Good 
agreement  between  calculated  and  observed  flutter  boundary  was 
obtained  for  most  all  cases.  The  flutter  frequencies  were 
found  in  between  the  support  pitch  frequency,  0.20  Hz  or  0.85 
Hz,  and  the  corresponding  wing  bending  frequencies.  The  best 
agreement  between  predicted  and  measured  frequencies  is  the 
model  with  the  win9*  For  roodel  with  the  [  3  0  2 /0  ]  s 
wing,  the  flutter  frequencies  were  somewhat  overpredicted  in 
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all  configurations,  while  those  of  the  [ — 1 52/0 ] s  wing  were 
somewhat  underpredicted.  A  possible  explanation  is  that  the 
overall  pitching  stiffness  was  not  correctly  predicted  in  the 
analytical  model  due  to  the  aeroelastic  twisting  and  untwisting 
of  the  bending-twi sting  coupled  composite  wings. 

This  well-correlated  analytical  model  was  attributed  to 
the  use  of  the  modified  strip  theory  and  the  measured 
aerodynamic  derivatives  as  described  in  Chapter  2.  The  use  of 
the  static  lift-curve-slope  distribution  in  wing  motion  proved 
to  be  a  good  approximation  because  of  the  low  flutter  reduced 
frequency,  k  -  0.02  to  0.05.  For  comparison  purposes,  the 
flutter  boundary  resulted  from  an  overall  finite  span 
correction  is  also  tabulated  in  Table  5.2.  The  overall 
correction  factor  is  given  by  Eq .  (5-223)  of  Ref.  35  as 

da  CL  "  AR  +  aQcosA/n  a0cosA  (5.1) 

where  ag  is  the  two  dimensional  lift-curve-slope  for 
incompressible  flow,  2n.  Results  show  that  Eq.  (5.1)  gives  a 
too  conservative  prediction. 

In  addition  to  the  flutter  boundary  and  flutter  frequency, 
the  measured  transient  aeroelastic  responses  were  also  compared 
with  the  predicted  responses.  Figures  5.3  to  5.6  show  the 
typical  predicted  and  measured  w-V  responses  on  the  modified 
support  system.  Results  of  the  other  tests  are  included  in 
Appendix  C.  In  general,  good  agreement  is  obtained  for  most 


Table  5.2  Comparison  of  Body  Preedoa  Flutter  Boundaries 
of  Two  Different  Finite  Span  Corrections 


Modified  Correction  by 
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10 
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cases.  Due  to  the  support  instability  inherent  in  the  original 
support  system,  the  transient  decay  measurements  incur  large 
scattering  in  wing  responses.  This  experimental  scattering 
could  arise  from  two  sources  when  support  instability  is 
present.  One  is  the  difficulty  in  experimentally  exciting  the 
wing  frequencies  in  this  region,  the  other  is  the  possible 
nonlinearity  in  flight  dynamics  (or  relaxed  wing  root  boundary 
condition).  Nevertheless,  results  from  the  original  and  the 
modified  support  system  converged  well  at  the  flutter  point. 
Because  of  the  better  controlled  testing  environment,  when  the 
plunging  mode  was  absent,  the  model  free  in  pitch  only 
exhibited  the  best  correlated  results. 

In  these  <*>-v  plots,  the  torsional-stall  and  bending-stall 
flutter  are  also  presented  without  analytical  correlation. 
This  stall  flutter  is  a  nonlinear  phenomenon  which  added  on  to 
the  body  freedom  flutter  predicted  by  the  present  linear 
analysis . 


5 . 2  Body  Freedom  Flutter  of  Completely  Unrestrained  Vehicle 

In  this  section,  the  aeroelastic  behavior  of  the 
completely  unrestrained  vehicle  will  be  compared  with  tunnel 
supported  configurations  analytically.  The  root  locus  analysis 
will  be  used  to  illustrate  the  transient  behavior  and  its 
instability  onset. 

In  a  root  locus  analysis,  the  complex  roots  of  the 
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characteristic  equation  are  plotted  on  the  complex  plane.  Each 
branch  of  complex  roots  starts  from  the  imaginary  (frequency) 
axis.  As  airspeed  increases,  the  complex  roots  move  on  the 
complex  plane  accordingly.  If  a  branch  of  complex  roots 
crosses  the  frequency  axis  and  moves  into  the  right  half  plane, 
then  instability  occurs.  Figures  5.7  to  5.10  show  root  locus 
plots  and  low  frequency  highlights  for  all  models  in  completely 
unrestrained  flight  condition.  In  these  figures,  each  symbol 
represents  a  3  m/sec  airspeed  increment. 

Comparison  of  these  root  locus  plots  shows  dramatic 
differences  in  transient  behavior  among  all  four  models.  For 
models  with  the  (02/90)s,  [302/0]s  and  l-152/0)s  win9s'  bodY 
freedom  flutter  instability  occurs,  since  the  pitch  mode 
becomes  unstable.  Examining  how  "fast",  (in  terms  of  speed 
increment),  the  pitch  mode  moves  into  the  right  half  plane,  the 
plots  show  that  the  model  with  [ — 1 5 2/° 1 s  wing  has  the  strongest 
instability,  while  the  model  with  [  3 0 2/0 ] s  wing  exhibits  a 
rather  weak  instability.  The  pitch  mode  of  the  [ 3  0  2/0  ]  s  winn 
is  eventually  stabilized  again  at  higher  speed  as  the  pitch 
branch  moves  back  to  the  left  half  plane.  The  model  with  the 
[152/0]s  wing  demonstrates  the  best  aeroelastically  tailored 
model,  since  the  pitch  mode  never  moves  into  the  right  half 
plane.  In  addition  to  the  body  freedom  flutter,  there  also 
exists  a  wing  bending- tor s ion  flutter  as  indicated  by  the 
crossing  of  the  torsional  mode  or  bending  mode  into  the  right 
half  plane.  It  should  be  noted  that  the  designations  of 
"pitch",  "bending",  and  "torsional"  modes  are  only  for 
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Fig.  5.7  Root  Locus  Plot,  Completely  Unrestrained  Model 
with  ( 0 2/90 ] s  Wing 
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Fig.  5.10  Root  Locus  Plot,  Completely  Unrestrained  Model 
with  (-15, /OJ  Wing  (cont'd) 
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convenience  notation.  As  airspeed  proceeds,  each  branch 
actually  comprises  contributions  of  various  degrees  of 
freedom.  For  example,  when  body  freedom  flutter  occurs,  the 
actual  flutter  mode  consists  of  model  pitch,  model  plunge,  wing 
bending,  and  wing  torsional  responses  as  indicated  in 
Table  4.4. 

Since  the  present  analytical  model  also  gives  the  true 
damping  (or  transient  decay  rate)  of  the  aeroelastic  system, 
the  influence  of  aeroelastic  tailoring  on  longitudinal  flying 
quality  can  be  assessed  in  a  preliminary  fashion  by  comparing 
the  damping  of  the  short  period  mode  (or  pitch  branch).  For 
example,  at  low  speed  range,  12  m/sec,  the  damping  (or  pitch 
branch  location)  is  relatively  unaffected  except  for  the  model 
with  the  [-152/0]s  wing.  Using  the  flutter  point  of  the  model 
with  102/90)s  win9*  18.5  m/s  ec,  as  the  design  point,  a 
comparison  between  the  tl52/0]g  and  [302/0]gwing  shows  small 
change  in  the  location  of  pitch  branch.  As  airspeed  increases, 
the  ( 302/0 ] s  wing  mildly  decreases  its  damping  while  the 
(152/0)s  wing  mildly  increases  its  damping.  At  higher 
airspeed,  the  model  with  (152/0]g  wing  actually  exhibits  higher 
damping  than  does  the  rigid  model. 

The  flutter  boundaries  and  frequencies  of  the  completely 
unrestrained  vehicle  are  also  summarized  in  Table  5.1. 
Comparisons  show  that  the  completely  unrestrained  vehicle  is 
closely  simulated  by  both  soft  (original)  and  stiff  (modified) 
support  systems.  Figures  5.11  to  5.14  compare  the  model 
transient  behavior  of  the  modified  support  system  with  that  of 
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the  completely  unrestrained  vehicle.  Again,  the  airspeed 
incremental  is  set  to  be  3  m/sec  for  the  root  locus  analyses. 
The  results  show  a  quick  convergence  to  the  behavior  of  the 
unrestrained  vehicle,  except  that  there  exists  a  support  plunge 
mode  in  the  tunnel  supported  configuration.  The  dynamic 
characteristics  of  these  support  plunge  modes  varies  with 
models.  For  models  with  the  {  0 2/9 0 } $  and  ( - 1 5 2/0 ] s  wings,  the 
support  plunge  mode  seems  to  be  more  active  than  is  the  model 
with  the  t 152/0 ] s  and  [302/0]s  wing.  In  general,  however,  the 
support  plunge  mode  tends  to  move  away  from  the  pitch  mode 
which  is  important  in  simulating  the  dynamics  of  a  completely 
unrestrained  vehicle. 

The  well-simulated  free  flight  condition  is  attributed  to 
the  carefully  designed  support  system.  The  support  stiffness, 
in  particular,  is  designed  to  be  well  separated  from  the 
suspected  range  of  flutter  frequencies,  and  support  frictions 
are  kept  to  a  minimum  to  prevent  contamination  of  the 
aerodynamic  damping  force  which  is  important  in  pinpointing  the 
body  freedom  flutter.  Additionally,  the  support  is  designed 
such  that  the  pitch  frequency  in  still  air  is  higher  than  the 
plunge  frequency,  in  order  to  avoid  the  tunnel  support  dynamic 
instability  reported  earlier. 

5 . 3  Effect  of  Some  Parameters  in  Body  Freedom  Flutter 


The  body  freedom  flutter  arises  from  the  interaction  of 
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the  aeroelastically  destiffened  wing  mode  with  the 
aerodynamically  stiffened  rigid  body  short  period 
(pitch-plunge)  mode.  Due  to  the  cantilever  wing  divergence 
tendency,  the  frequency  of  the  wing  mode  decreases  with 
airspeed,  while  the  short  period  mode  frequency  increases  with 
airspeed.  As  a  result,  there  exists  a  dynamic  instability 
characterized  by  the  phenomenon  of  frequency  coalescence.  The 
exact  occurrence  of  body  freedom  flutter  will  be  functions  of 
those  parameters  pertinent  to  the  cantilever  wing  divergence 
and  aircraft  flight  dynamics. 

The  so-called  "aeroelastic  destif fening" ,  an  underlying 
mechanism  of  the  aeroelastic  wing  divergence,  is  not  only  a 
function  of  the  wing  frequencies  (or  stiffness)  but  also  a 
function  of  the  elastic  coupling  for  an  aeroelastically 
tailored  wing.  Effect  of  the  elastic  coupling  has  been 
summarized  in  Sec.  5.1. a.  To  isolate  these  two  effects,  the 


influence  of  wing  frequencies  will  be  discussed  on  the  model 
with  the  [02/90]s  wing.  On  the  other  hand,  the  flight  dynamics 
effect  in  symmetric  flight  is  primarily  the  short  period  mode 
effect.  The  important  parameters  involved  in  flight  dynamics 
are  the  "rigid"  aircraft  static  margin  and  inertia  coupling  of 
wing  elastic  modes  and  rigid  body  modes. 

With  this  verified  analytical  model,  the  aeroelastic 
behavior  of  forward  swept  wing  aircraft  can  be  further 
understood  from  discussions  of  the  aforementioned  parameters. 
It  should  be  noted  that,  for  the  purposes  of  parametric 
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However,  the  intended  parametric  configuration  may  not  be 
physically  possible.  In  other  words,  a  real  design  parameter, 
e.g.,  wing  root  position,  may  affect  several  parameters,  i.e., 
aircraft  static  margin  and  inertia  coupling. 


( a )  Influence  of  wing  frequencies 


As  mentioned  earlier,  the  influence  of  wing  frequencies 
will  be  discussed  only  for  the  model  with  the  [02/90]s  win9- 
Using  the  current  analytical  model,  the  parametric  predictions 
are  illustrated  in  Figure  5.15.  The  wing  fundamental  bending 
frequency  is  changed  by  artificially  varying  the  value  of 
the  AS4/3501-6  graphite/epoxy  while  the  other  material 
properties  remain  unchanged.  Similarly,  the  wing  torsional 
frequency  is  changed  by  parametrically  varying  the  GLT  value 
only. 


The  resulting  body  freedom  flutter  speeds  show  the 
decreasing  trend  with  lowered  wing  bending  frequency.  The 
outcome  is  expected  as  the  possibility  of  frequency  coalescence 
and  wing  destiffening  are  both  enhanced  by  lowered  wing 
fundamental  bending  frequency.  On  the  other  hand,  an  important 
but  less  familiar  finding  shows  that  the  body  freedom  flutter 
speed  also  decreases  with  the  lowered  wing  torsional  frequency 
which  is  presumed  to  be  less  relevant  to  the  interaction 
mechanism.  For  a  cantilever  wing,  the  torsional  stiffness 
(torsional  frequency)  is  known  to  be  a  deciding  factor  in 
torsional  divergence.  This  result  reconfirms  an  important 
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mechanism  in  the  body  freedom  flutter,  i.e.,  the  interaction  of 
the  destiffened  wing  mode  (or  wing  "divergence"  mode).  This 
parametric  study,  identifies  an  ancestry  root  of  wing 
"divergence"  mode  in  the  body  freedom  flutter  instability. 
This  result  also  suggests  that  a  realistic  body  freedom  flutter 
analysis  should  at  least  include  both  wing  bending  and 
torsional  freedoms  to  prevent  the  erroneous  conclusions. 

( b )  Effect  of  "rigid"  aircraft  static  margin 

The  presence  of  flight  dynamics  interaction  can  be 

physically  seen  as  the  modification  of  boundary  conditions  of 

the  destiffened  wing  mode.  These  modified  boundary  conditions 

include  both  "stiffness"  and  "damping”  effects  resulting  from 

the  aircraft  aerodynamic  characteristics.  The  "rigid"  aircraft 

static  margin,  determined  by  the  relative  positions  of  aircraft 

center  of  mass  and  lifting  surfaces,  is  conveniently  used  to 

characterize  these  aerodynamic  properties.  in  the  aeroelastic 

analysis,  however,  varied  static  margin  accomplished  by  lifting 

surfaces  relocation  will  also  affect  the  degree  of  inertia 

coupling  between  the  elastic  and  rigid  body  motions.  For  the 

present  analytical  model,  the  use  of  rigid  model  aerodynamic 

derivatives  makes  an  isolated  parametric  study  of  static  margin 

effect  possible.  In  other  words,  the  static  margin,  ratio  of 

-CM  to  ,  will  be  varied  by  changing  the  CM  accordingly 
a  a  a 

while  maintaining  the  value  the  same. 

a 

This  parametric  study  is  equivalent  to  relocating  canard 
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surface  only,  but  neglecting  its  inertia  effect  which  is  true 
for  a  rigid  canard  model.  To  hypothesize  a  wide  range  of 
static  margin,  0%  to  60%  of  streamwise  chord  length,  the  actual 
canard  position  may  be  physically  off  the  fuselage  or 
interfering  with  other  lifting  surfaces,  but  these  effects  are 
ignored  here. 

Figure  5.16  illustrates  the  static  margin  effect  on  body 
freedom  flutter.  Both  flutter  speed  and  frequency  are 
increased  as  static  margin  increases.  This  effect  is 
consistent  for  all  wings  but  with  varied  degrees  of  sensitivity 
to  static  margin.  For  example,  when  static  margin  is  increased 
from  0%  to  60%,  the  model  with  [0-/90]  wing  has  only  2% 
increase  in  flutter  speed,  while  the  [ 3 0 2/0 1 s  wing  registered  a 
13%  increase.  To  demonstrate  the  possibility  of  dramatic 
influence  of  static  margin,  the  model  with  [102/0]  wing  is 
included  in  this  study.  The  results  show  that  the  body  freedom 
flutter  can  be  significantly  delayed  or  eliminated  by  the 
increased  static  margin.  At  the  first  glance,  the  results  seem 
to  contrast  to  the  physical  interpretation  of  the  adverse 
effect  of  "aerodynamic  stiffening".  Since  the  aerodynamic 
stiffness  is  proportional  to  the  aircraft  static  margin,  or 

,  the  increase  of  static  margin  should  promote  the 
a 

possibility  of  earlier  frequency  coalescence,  and  result  in  a 
lowered  flutter  speed.  However,  the  increased  static  margin 
will  also  simultaneously  increase  the  coupled  aerodynamic 
damping  and  its  stabilizing  effect.  As  a  result,  the 
aerodynamic  damping  is  significant  enough  to  delay  the 


STATIC  MARGIN  ( t  of  streamwise  chord) 


It  is 


destabilising  trend  of  frequency  coalescence.  It  is 
interesting  to  note  that  the  analysis  with  mechanical  pitching 
stiffness  can  demonstrate  the  isolated  pitching  stiffness 
effect.  As  expected,  the  more  the  added  mechanical  pitching 
stiffness,  the  lower  the  body  freedom  flutter  speed. 

For  aircraft  with  negative  static  margin,  the  aircraft 
static  instability  is  more  critical  than  is  the  body  freedom 
flutter.  The  investigation  of  body  freedom  flutter  then 
becomes  meaningless  unless  a  stability  augmentation  system 
(SAS)  is  implemented  and  incorpoated  in  the  analytical  model. 
In  Ref.  24,  body  freedom  flutter  was  found  to  be  dependent  on 
the  closed  loop  characteristics  of  the  SAS,  and  similar  to  the 
behavior  of  a  statically  stable  aircraft. 
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(c)  Inertia  Coupling  Effect 

In  addition  to  the  "stiffness"  and  "damping"  effects  S 
introduced  by  aircraft  flight  dynamics,  the  presence  of  rigid 
body  freedoms  also  provides  the  inertia  coupling  effect  (or 
inertia  relief)  between  wing  elastic  modes  and  rigid  body 
modes.  With  the  present  analytical  model,  this  coupling  effect 

.v 

is  examined  by  varying  wing  root  position  while  keeping  the 
"rigid"  aircraft  static  margin  the  same.  Figure  5.17  shows 

Cs 

that  body  freedom  flutter  decreases  as  the  wing  moves  aft.  °  , 
This  is  due  to  the  increased  pitch  inertia  coupling,  nu,  in  ^ 
present  study,  as  evident  by  the  corresponding  vibration  mode. 

Another  straightforward  indication  of  inertia  coupling  is  ■*'  ! 


I’.  HOOT  POSITION  AFT  AIRCRAFT  C.G.  U  ot  lusel.t,.  length) 
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the  wing  to  aircraft  mass  ratio.  Practically,  any  change  in 
wing  mass  ratio  will  affect  both  pitch  and  plunge  inertia 
coupling  simutaneously .  For  the  purposes  of  isolated 
parametric  study,  the  wing  mass  and  pitch  inertia  ratios  will 
be  varied  seperately  by  changing  the  overall  aircraft  mass  and 
pitch  inertia  respectively  without  changing  wing  mass. 
Figure  5.18  shows  the  adverse  effect  of  increased  wing  pitch 
inertia  ratio  (decreased  aircraft  pitch  inertia),  or  increased 
pitch  inertia  coupling,  m2  3//m2  2  ’  5.19  shows  the 
favorable  effect  of  increased  wing  mass  ratio  (decreased 
aircraft  mass),  or  increased  plunge  inertia  coupling,  m^./m^. 
For  the  present  wind  tunnel  model,  the  measured  wing  mass  and 
pitch  inertia  ratios  are  5.6%  and  3.11  respectively. 


4 
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I 
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5 . 4  Aeroelastic  Effect  on  Support  Instability 


The  testing  and  analytical  results  have  shown  that  the 
model  support  instability  interferes  very  little  with  the  body 
freedom  flutter.  It  is  also  of  interest  to  further  explore  the 
mechanism  involved,  and  to  foster  the  understanding  of  both 
instabilities. 

From  the  analysis,  it  is  shown  that  the  lower  limit  of  the 
support  instability  zone  is  dictated  by  the  support  pitch 
damping  force,  while  the  upper  limit  is  determined  by  the 
support  plunge  damping  force,  when  the  wing  flexibility  is 
taken  into  account,  there  are  additional  aerodynamic  damping 
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Fig.  5.19  Plunge  Inertia  Coupling  Effect  on 
Body  Freedom  Flutter 


contributions  from  the  wing  elastic  motions.  In  general,  the 
wing  motions  can  be  decomposed  into  wing  bending  and  twisting 
components.  The  aerodynamic  damping  incurred  by  the  wing 
bending  motion  is  usually  significant  and  additive  to  the 
overall  plunge  damping.  The  wing  bending  contribution  to  the 
overall  pitch  damping  is  less  significant  due  to  the  proximity 
in  location  of  pitching  axle  and  wing  tip  which  would  have  the 
most  wing  bending  motion.  On  the  other  hand,  the  wing  twisting 
motion  is  relatively  small  in  amplitude,  and  it's  contribution 
to  the  overall  plunge  and  pitch  damping  about  model  pitching 
axle  is  negligible.  As  a  result,  the  aeroelastic  effect  on 
support  instabiity  come  mainly  from  the  wing  bending  motion  and 
its  contribution  to  the  overall  plunge  damping  force.  This 
effect  was  observed  in  wind  tunnel  tests  and  predicted  from  the 
analytical  model.  The  results  are  illustrated  in  Figure  5.20. 

A  close  examination  of  Figure  5.20  shows  that  the  upper 
instability  boundary,  31  ra/sec  for  "rigid"  model,  primarily 
follows  the  trend  of  the  body  freedom  flutter.  This  is  because 
the  body  freedom  flutter  mode  always  contains  significant  wing 
bending  motion.  Without  the  presence  of  any  aeroelastic 
instability,  the  aeroelastic  response  is  usually  small.  When 
the  body  freedom  flutter  instability  is  approached,  the  wing 
bending  response  becomes  significant  and  the  overall  plunge 
damping  is  eventually  dominated  by  the  wing  bending  motion. 
The  upper  boundary  of  support  instability  is  therefore  dictated 
by  the  trend  of  body  freedom  flutter  boundary.  in  general,  any 
aeroelastic  instability  involving  wing  bending  motion  such  as 
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bending-torsion  wing  flutter  will  have  the  same  effect  as  body 
freedom  flutter  does.  The  support  instability  of  the  model 
with  the  ( 1 5 2/0 ) s  wing,  which  is  "optimum"  in  term  of  body 
freedom  flutter,  is  actually  controlled  by  the  critical 
bending-torsion  flutter  instead  of  the  body  freedom  flutter. 

Since  the  body  freedom  flutter  boundary  is  well  predicted 
in  the  present  analysis,  the  upper  limit  of  support  instabi  lity 
is  subsequently  shown  to  be  in  good  agreement  with  the  testing 
results  except  for  the  model  with  the  [  3 0 2/0  ] g  wing.  in  this 
case,  the  support  instability  terminates  at  23  m/sec  compared 
with  the  predicted  value  of  16  ra/sec.  A  close  examination  of 
the  testing  records  shows  that,  between  the  upper  limit  of 
support  instability  and  the  body  freedom  flutter  onset,  the 
system  was  so  lightly  damped  that  any  small  disturbance  would 
cause  the  model  to  plunge  significantly  at  such  airspeed. 
During  wind  tunnel  tests,  it  was  difficult  to  distinguish  such 
responses  from  the  support  instability  even  though  it  may  be 
due  to  a  turbulent  gust. 

Unlike  the  upper  limit,  the  lower  limit  of  support 
instability  is  governed  by  the  support  pitch  damping  force. 
This  involves  very  little  aeroelastic  effect,  and  the  accuracy 
of  the  lower  limit  depends  on  the  correct  estimation  of  the 
support  pitch  damping.  Because  the  model  was  heavily  damped  in 
pitch  freedom  in  the  original  support  system,  only  about  one 
and  one  half  cycle  of  a  step  response  was  observable,  and  this 
led  to  some  scattering  and  uncertainty  in  determining  the 
damping  ratio.  Nevertheless  the  slight  underpcedictions  in 


lowec  Unit  of  support  instability  are  consistent  among  all 
wings  and  within  the  experimental  error  range. 

Since  the  body  freedom  flutter  has  a  "healing"  effect  on 
the  support  instability,  which  occurs  first,  the  latter  is  seen 
to  have  little  opportunity  to  affect  the  former.  This 
conclusion  of  a  "one  way"  effect  is  proved  from  the  analysis 
and  testing  results. 
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CHAPTER  SIX 

CONCLUSIONS  AND  RECOMMENDATIONS 

An  analytical  and  experimental  investigation  was  made  of 
the  aeroelastic  flutter  and  divergence  behavior  of  a  forward 
swept  graphite/epoxy  wing  aircraft  with  rigid  body  pitch  and 
plunge  freedoms.  Experimental  wind  tunnel  data  were  obtained 
for  a  generic,  full-span,  30®  forward  swept  wing,  18%  static 
margin  aircraft  model  in  a  low  speed  wind  tunnel.  The  free 
flying  behavior  was  simulated  on  a  softly  restrained  vertical 
rod  support  system.  Based  on  the  results  of  this 
investigation,  the  following  conclusions  are  drawn: 

1.  The  tests  revealed  large  variations  in  aeroelastic 

behavior  for  the  aircraft  model  with  four  different  ply 
layup  wings.  Body  freedom  flutter  and  bending-torsion 
wing  flutter  were  observed  on  the  free  flying  model,  and 
regions  of  their  occurrence  were  mapped  out. 

2.  The  forward  swept  wing  aircraft  model  developed  body 
freedom  flutter  rather  than  divergence  when  rigid  body 
freedoms  were  present.  This  body  freedom  flutter  occurred 
at  a  lower  airspeed  than  cantilever  wing  divergence,  and 
arose  from  the  interaction  of  the  destiffened  wing  mode 
and  the  aircraft  short  period  mode.  The  presence  of 
flight  dynamics  effects  modified  the  boundary  conditions 
of  the  destiffened  wing  mode,  and  ultimately  changed  the 
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instability  nature  from  static  divergence  to  dynamic 
flutter. 


The  aeroelastic  tailoring  effect  on  body  freedom  flutter 
was  found  to  follow  the  trend  of  cantilever  wing 
divergence  tailoring.  The  model  with  the  wing 
demonstrated  the  best  tailored  aeroelastic  behavior,  and 
bending- tor sion  wing  flutter  was  encountered  there  instead 
of  body  freedom  flutter.  The  presence  of  flight  dynamics 
was  found  to  have  little  effect  on  this  bending-torsion 
type  of  wing  flutter. 


A  support  related  dynamic  instability,  involving  the 
model's  pitch  and  plunge  motions,  was  also  encountered  at  ^ 
low  airspeeds.  This  support  instability  could  be 
eliminated  by  making  the  still  air  model  pitch  frequency  [-1 

r  % 

greater  than  the  model  plunge  frequency,  thereby  avoiding 
a  possible  frequency  coalescence  at  certain  airspeed.  It 

*  * 

was  also  found  that  this  support  instability  was  banded  in 
airspeed,  and  caused  almost  no  effect  on  the  body  freedom  >• 
flutter  boundaries,  if  the  model  could  fly  through  the 
instability  zone. 


For  the  model  free  in  pitch  only,  body  freedom  flutter  was 
again  encountered,  but  at  a  lower  airspeed,  and  lower 
flutter  frequency.  The  presence  of  the  plunge  freedom 
apparently  played  an  important  and  favorable  role  in  body 
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freedom  flutter,  and  should  not  be  overlooked  in  body 
freedom  flutter  analysis  or  test. 

6.  Tests  conducted  on  the  cantilever  wings,  i.e.,  the 

aircraft  model  locked  in  plunge  motion  and  restrained  in 
pitch  motion  by  the  snubber  cables,  revealed  the  same 
aeroelastic  behavior  as  noted  earlier  in  Ref.  10.  Linear 
wing  divergence  and  bending-torsion  wing  flutter  occured 
at  low  angles  of  attack,  and  nonlinear  torsional-stall 
flutter  and  bending-stall  flutter  occured  with  lowered 
airspeed  at  the  higher  angles  of  attack.  Depending  on  the 
wing  layups,  the  static  wing  divergence  condition  could 
change  to  either  a  torsional-stall  flutter  or 

bending-stall  flutter  as  wing  root  angle  of  attack  was 

increased.  The  flutter  prone  wing  always  exhibited 
torsional-stall  flutter  at  the  higher  angles  of  attack. 

7.  An  accurate  and  efficient  linear  flutter  analysis  was 

developed  by  using  Rayleigh-Ri tz  approximation,  modified 
strip  theory  in  transient  representation,  and  measured 
aircraft  aerodynamic  derivatives.  Good  agreement  was 

obtained  for  all  cases  in  both  the  flutter  boundaries  and 

transient  behavior.  The  use  of  static  lift-curve-slope 
distribution  in  modified  strip  theory  was  shown  to  be  a 
good  approximation  because  of  the  low  body  freedom  flutter 
reduced  frequency,  k«  0.02  to  0.05. 
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8.  parametric  studies  based  on  the  analysis  developed  here 
showed  that  both  the  wing  bending  and  torsional 
frequencies  were  critical  to  the  body  freedom  flutter  as 
well  as  to  the  cantilever  wing  divergence.  The  increased 
rigid  aircraft  static  margin  showed  a  favorable  effect  in 
body  freedom  flutter  speed  for  all  wings  analyzed.  The 
plunge  inertia  coupling  of  the  wing  elastic  modes  and 
rigid  body  modes  shows  a  favorable  effect,  whereas  the 
pitch  inertia  coupling  shows  an  adverse  result. 


3 


The  present  study  has  provided  analytical  and  experimental 


insight  into  the  actual  aeroelastic  behavior  of  a  forward  swept 
wing  aircraft  in  free  flight.  However,  there  are  several  other 


important  analytical  and  experimental  aspects  related  to  this 


subject  that  require  additional  investigation.  Recommendations 


for  further  study  are  given  below. 


From  an  analytical  viewpoint, 


effects  of 


compressibility  should  be  evaluated  to  complete  the 


parametric  trend  prediction, 


Regarding  aerodynamics,  the  effect  of  nonlinearity  at 
higher  angle  of  attack,  which  cause  a  transition  from 


linear  behavior  to  torsional-stall  and  bending-stall 


flutter,  deserves  further  investigation. 


Based  on  this  well-verified  aeroelastic  plant  model,  there 
is  a  desire  to  demonstrate,  through  analysis  and  wind 
tunnel  tests,  the  integration  of  aeroelastic  tailoring 
technique  and  active  control  technology  for  instability 
suppression.  The  inclusion  of  rigid  body  freedom  are 
certainly  desired  here. 

Analyses  and  wind  tunnel  tests  also  need  to  be  conducted 
to  evaluate  the  antisymmetric  aeroelastic  behavior  of  a 
forward  swept  wing  aircraft  in  free  flight.  From  there, 
the  study  can  be  further  extended  to  the  subject  of 
oblique  wing  design. 
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APPENDIX  A 

STRUCTURAL  INERTIA  AND  STIFFNESS  COEFFICIENTS 


A . 1  INERTIA  COEFFICIENTS 
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M  is  the  wing  mass,  and 
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constants  as  tabulated  in  Table  A.l.  It  should  be  noted  that 
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the  chosen  mode  shapes  uncouple  the  elastic  degrees  of  freedom, 
but  there  still  exists  inertia  coupling  between  the  rigid  body 
and  elastic  modes,  i.e.,  ra, ^ »  nu  ^  etc. 


A. 2  STIFFNESS  COEFFICIENTS 


Similarly,  the  stiffness  coefficients  given  by  Eq.  (2.8 
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It  should  be  noted  that  the  torsional  stiffness  terms,  k44  and 
kgg,  are  modified  for  the  additional  warping  stiffness  by 
correction  factors,  k^t  and  k2t‘  is  because  that  the 
assumed  St.  Venant's  torsional  modes  are  free  to  warp  which  is 
incompatible  with  the  root  boundary  condition  of  the  plate 
model.  Reference  43  detailed  these  warping  corrections  based 
on  the  partial  Ritz  solution  of  a  symmetric  laminated  plate 
without  bending-twisting  coupling,  i.e.,  D^g-0  and  D26"^‘ 
Figure  A.l,  reproduced  from  Ref.  10,  shows  k^t  and  k£t  versus 
the  effective  aspect  ratio,  0,  which  is  defined  as 
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Table  A.l  Mode  Shape  integrals  (cont'd) 


u-  l2  J  (  ♦g  *4)  dx  -  -6.700 
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APPENDIX  B 

TRANSIENT  AERODYNAMIC  COEFFICIENTS 


V  , 


The  two-dimensional  incompressible  harmonic  airloads, 
including  the  effect  of  camber  oscillation,  was  given  by 
Spielberg  in  Ref.  31.  Using  the  sign  conventions  and  cambering 
mode  (chordwise  bending  mode)  defined  in  Chapter  2,  this 
result,  about  mid-chord  location,  was  rederived  and  are  given 
in  the  following  expressions  as 
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To  obtain  the  transient  airload  expressions,  a  key  step  is 
to  generalize  the  airload  expressions  from  the  frequency  axis 
into  the  complex  variable  domain  [45,46].  For  the  aerodynamic 
model  used  in  the  presnet  study,  one  can  first  substitutes 
nondimensional  Laplace  variable,  p,  for  jk  and  generalized 
Theodorsen  function,  C(p),  for  C ( k )  in  Eq.  (B.l),  i.e., 


L-  2>lpV2bcosA([A2ap2>Aup+A0a+(A4ap+A3a)C(p)]  g 


[A2bP2+AlbP+A0b+ ( A4bP+A3b  > C ( P  > ]  a 
[A2cP2+AlcP+A0c"(A4cP+A3c)C(P)]  b  } 

M-  InoV2b2COSA|[A2dp2i.Aldp>A0d^(A4dp+A3d)C(p)]  £ 


[A2.P: 


+AleP>A0eMA4eP^A3e)C(p)J  * 


(  B  .  2  ) 


[A2fP  *A1  fP*A0  f  *  (  A4  f  P+A3  f  J  C  (  P  * .  b 


-185- 


N-  2R0V2bcO5A|[A2gp2+Algp+A0g+(A4gp+A3g)C(p)]  | 

[A2hp2*AlhP-A0h+ ( A4hPtA3h 1 C  < P » ]  * 
[A2iP2+AliP*A0i-< A4tP-A3i >C(P>]  I  } 

where 


4a" 

-1 

A3a“ 

0 

A2a“ 

i 

'7 

ai»- 

0 

> 

o 

0) 

• 

o 

4b“ 

1 

7 

A3b“ 

1 

A2b" 

0 

Aib- 

1 

7 

A0b"  0 

4c" 

1 

T 

A3c" 

-1 

A2c" 

1 

77 

Alc" 

0 

A0c"  0 

4d* 

0 

A3d“ 

1 

* 7 

A2d* 

0 

Ald- 

0 

A0d"  0 

4e“ 

l 

4 

A3e" 

l 

7 

A2e" 

1 

“17 

Ale- 

l 

“7 

A0e"  0 

4  f  “ 

1 

'17 

A3  f  " 

1 

'7 

A2  f  " 

0 

Alf- 

l 

4 

A0f"  7 

4g“ 

l 

A3g" 

0 

A2g* 

i 

27 

A19- 

0 

A0g"  0 

4h“ 

l 

17 

A3h" 

l 

Z 

A2h" 

0 

Alh’ 

1 

~z 

A0h"  0 

4  i" 

l 

'T& 

A3i" 

l 

A2  i" 

l 

"77 

AU* 

0 

A0i”  7 

The  next  step  is  to  approximate  these  motion-dependent  airloads 
by  functions  of  displacements,  their  derivatives,  and  finite 
number  of  lag  terms  as 
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The  coefficients  of  Eq.  (B.3),  e.g.,  B2a>  bia»  bqa'  and  BnA' 
can  be  thought  of  as  the  aerodynamic  inertia,  damping, 
stiffness,  and  lag  terms. 

To  determine  those  coefficients  of  Eq.  (B.3),  one  can 
usually  fits  Eq.  (B.2)  to  a  table  of  harmonic  airloads  computed 
for  a  chosen  set  of  reduced  frequencies.  In  the  case  of 
two-dimensional  incompressible  flow,  this  curve  fitting  can  be 
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actually  performed  for  the  generalized  Theodorsen  function, 
C(p) ,  instead  of  the  entired  airloads  calculated  from 
Eq.  (B.2). 

For  numerical  fitting,  the  number  of  lag  terms  and 
corresponding  pole  locations  are  somewhat  arbitrary,  but 
influence  the  accuracy  of  the  functional  fit.  In  many  cases, 
the  least  squares  error  function  minimization  was  often  used  in 
finding  these  coefficients  or  pole  locations  [47,48]. 
Typically,  numerical  experimentation  is  used  to  obtain 
reasonable  fits  in  the  anticipated  reduced  frequency  range  with 
as  few  lag  terms  as  possible. 

For  the  present  investigation,  a  single-lag  approximation 
was  chosen  with  pole  located  at  -0.15,  i.e., 
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where  b^-  0.55  and  bg-  0.15.  Placing  Eq.  (B.4)  into  (B.2),  and 
rearranging  terms  according  to  Eq.  (B.3)  with  L«3  results  in 
the  following  set  of  coefficients 
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Introducing  new  augmented  state  variables,  y^,  such  that 
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and  placing  these  augmented  state  variables  into  Eq .  (B.5) 


leads  to  the  following  mixed  variable  expressions 


-13  j- 


*  2 
L-  2«pV  bcosA 

{[s2«p2*»i.p*Bo.]£1\-Tr  *  *3.^  \^r 

r  2  i  7  7 

LB2bp  *Blbp*B0bJ  X  9i  “H  *  B  3b  X  ®i  yi 

(  B  .  7  ) 

[B2cp2'Blcp*80c]  X^iTT  *  S3cX1  ‘iTT  j 

2ixpV^b2cosA  •••  etc 
N-  2npV^bcosA  ••••  etc 

Taking  the  inverse  Laplace  transform  of  Eq .  (B.7)  and  (B.6) 
results  in  a  two-dimensional,  incompressible  airload 
expressions  in  time  domain  and  governing  equations  for 
aerodynamic  lag  terms.  To  approximate  the  aerodynamic  effect 
of  a  finite  span  wing,  these  sectional  airload  expressions  will 
be  modified  by  substituting  the  static  1 i f t-cur ve-s 1  ope 
distribution,  ,  for  2ncosA  such  that 
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